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Abstract 



The generahy accepted phase diagrams for the discrete Z^r spin 
niodels in two dimensions imply the existence of certain renormali- 
sation group flows, both between conformal field theories and into a 
massive phase. Integral equations are proposed to describe these flows, 
and some properties of their solutions are discussed. The infrared be- 
haviour in massless and massive directions is analysed in detail, and 
the techniques used are applied to a number of other models. 
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1 Introduction 



A planar vector model in which an 0(2) symmetry has been broken down 
to a subgroup can be defined by the reduced Hamiltonian 



The pairs {ij) are nearest-neighbour sites of the two-dimensional lattice on 
which the angular spins 0, live; these interact via a 27r-periodic, even function 
V{9). The strength of the explicit symmetry-breaking is governed by the 
coupling h. 

When h=0, continuous symmetry is regained and this forbids the forma- 
tion of an ordered phase, no matter how low the temperature T becomes. 
Instead, the model undergoes a Kosterlitz-Thouless transition as T decreases 
through some critical value Tc. Below this temperature, the completely dis- 
ordered high-temperature phase, in which the typical spin configuration is 
dominated by vortices, is replaced by a massless phase in which the vortices 
are suppressed and correlations decay algebraically. For h 7^ 0, this picture 
must change: since residual symmetry is discrete there is no obstruction 
to its spontaneous breaking at low enough temperatures, and a rigorous 
result of Frohlich and Lieb asserts that this does indeed happen Q. De- 
spite this fact, for N > 5 the massless phase is not completely lost: Jose et 
al Q showed that the perturbation is irrelevant above a temperature T^' 
which, for > 5, lies below Tc- A small non-zero value of h cannot change 
the h=0 critical behaviour until T falls below T^', and so the Kosterlitz- 
Thouless phase survives as an intermediate stage between massive high- 
and low- temperature regimes. 

Another situation amenable to analysis is the limit h ^ cxd: this serves 
to pin the spins to N discrete values, and results in a 'pure' model. Uni- 
versality would suggest the continuing validity of the picture just outlined, 
and this expectation was confirmed by Elitzur et al in 1979 |^. However, 
both this and the earlier work of Jose et al concentrated on the Villain 
form IQ] of the interaction V{6), and it is important to know how much 
influence this choice has on the nature of the phases and transitions ex- 
hibited by the model. For the pure case, it is feasible to investigate 
this question directly, since the full phase space is finite ([A/2]) dimen- 
sional: a particular system is specified by the N reduced energy differences 
Vr = V{2Trr/N)—V{0), subject to the condition Vr = VN-r- The general 
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features of the resulting phase diagrams have been the subject of much at- 
tention since the early 1980's, motivated both by the intrinsic interest of the 
problem and by the analogies between such two-dimensional spin systems 
and certain four-dimensional gauge theories 

For each N, a duality transformation can be defined as follows Q. First 
parametrise the couplings as Xr = exp{—Vr). Then the summation over the 
spins 9i can be replaced by one over dual spins 9^ defined on the sites i of the 
dual lattice, so long as the set of couplings {xr} is simultaneously replaced 
by the dual set {xr}'- 

The effect is to exchange the low and high temperature fixed points {xr=0 
and Xr=l respectively), while leaving an [A^/4] -dimensional hyperplane in- 
variant - systems on this hyperplane are said to be self-dual. 

Setting a^; = e*'^^' and /i[ = e*'^' , the various phases can be characterised 
by the order parameters (a^) and their duals, the 'disorder parameters' 
{jj,'-). In a completely ordered phase, (cr'^) 7^ 0, (^') = 0; in a completely 
disordered phase (cr^) = 0, (/i') 7^ 0; while both vanish in a Kosterlitz- 
Thouless phase. In addition, if is not prime there are 'partially ordered' 
phases where expectation values of certain powers (a^) or (fi^) are nonzero 
even though others vanish Q. However such phases, and the corresponding 
partial-ordering transitions, are found away from the region that will be of 
interest below. 

The way that the three phases which are important for this paper fit 
together can be illustrated with the case N=5, shown in figure |l|. The 
transition from the ordered to the disordered phase can happen in one of 
two ways: either via an intermediate massless phase, as along the Villain line 
(6), or else directly through a single first-order transition, as exemplified by 
the behaviour along the Potts line (a). Clearly, the points C and C on the 
figure are rather special: they mark the bifurcations of the line of first-order 
transitions into pairs of Kosterlitz-Thouless transitions, and the beginnings 
of the massless regions. These points are thought to be special for other 
reasons too. In 1982 Fateev and Zamolodchikov identified particular sets of 
self-dual Boltzman weights for the general model which satisfy the star- 
triangle relations, and conjectured that they were precisely the critical points 
C and C for N=5, or their generalisations to higher |^]. Subsequently, in 
1985, they found for each N a ^Tv-symmetric conformal field theory that was 
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a natural candidate for the field-theory limit of these same special points . 




Figure 1: Phases for the general ^5 model: the Kosterlitz-Thouless regions 
1 and 2 are massless, while regions 3 (ordered) and 4 (completely disordered) 
are massive. The line AB is self-dual; a labels the line of Potts models, and 
h is the Villain line. 



To verify these claims, two distinct questions must be addressed. First, 
the criticality of the Fateev-Zamolodchikov Boltzman weights should be es- 
tablished and their critical exponents matched with those of the correspond- 
ing ^AT-symmetric conformal field theory: this has now been done, in both 



numerical and analytical |11] studies, with complete agreement with ex- 
pectations. However, on its own this information is not enough to locate the 
models relative to the opening of the first-order transition into a massless 
phase. To answer this second question, neighbouring regions of the phase 
diagram must be examined. Alcaraz |jl^ found numerical results consistent 
with the conjecture for the ^5 case, but beyond this, direct evidence from 
the lattice model is hard to come by. But the question can also be studied 
in the continuum, using the techniques pioneered by A.B.Zamolodchikov in 
his work on perturbed conformal field theories If the system described 
by the ^Ar-symmetric conformal field theory does indeed lie on the bor- 
der of a Kosterlitz-Thouless phase, then a suitable perturbation will shift 
it either into a massive phase, or else into the massless Kosterlitz-Thouless 
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region. Correspondingly, the conformal theory should admit perturbations 
'S'cFT — *■ 'S'cFT+A / e d?x which provoke both massive and massless renormal- 
isation group flows, according to the sign of the coupling A. The perturbing 
operator e should be ^Tv-symmetric, since the initial shift of the Hamiltonian 
is within the phase diagram of the 2jv models; and taking e to be self-dual 
will ensure that any massive direction of the perturbation is onto the surface 
of first-order transitions. The massless direction, shifting the model into a 
Kosterlitz-Thouless phase, should lead to a flow from c = 2(N—l)/(N+2), 
the central charge of the ^Tv-symmetric conformal field theory, to c = 1, the 
central charge of the various one-component Gaussian models. 

If e can be chosen so that the perturbed theory is integrable, then past 
experience of a method known as the Thermodynamic Bethe Ansatz (TBA) 
suggests that there should exist a set of integral equations encoding the 
exact flow of the ground-state energy E{R) of the system on a cylinder of 
circumference -R, from i?=0 all the way through to R=oo. Equations of 
this type were initially derived for theories for which a conjectured exact 
S-matrix was already available [13, 14 1, but it has since proved possible to 
make educated guesses for equations to describe many ground-state energy 
flows even in the absence of direct derivations |15, 16|. 

The purpose of this paper is to propose just such sets of equations, per- 
tinent to self-dual flows from the ^^v-symmetric conformal field theory into 
either a massive or a Kosterlitz-Thouless phase. The systems themselves are 
given in the next section, while section 3 shows that at small R they match 
a suitable perturbation already shown by Fateev to be integrable [0]. The 
discussion then turns to the large R behaviour, starting with the massless 
direction. In section 4.1, this is shown to be compatible with a self-dual 
flow into the Kosterlitz-Thouless region. We were able to extract some ex- 
act infrared information directly from the TBA, and this is explained in 
section 4.2. The technique is applied in a number of contexts, with results 
that reflect the connections between the various models associated with c=l 
conformal field theories. In the massive direction, described in section 4.3, 
a detailed examination of the asymptotics allows us to conjecture a pattern 
of massive kinks linking A'^-l-1 degenerate vacua, as expected on a surface 
of first-order transitions where phase coexistence occurs. The concluding 
section 5 outlines some questions that might merit further study. 

A number of the more technical details have been relegated to a series of 
appendices: the integral kernels, Y-systems and dilogarithm sum rules used 
in the main text can be found in appendices A, B and C, while numerical 
results for the ultraviolet asymptotics are reported in appendix D. 
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2 TBA equations 



In this section, the equations will be introduced, leaving a detailed study of 
their properties until later. Some are already in the literature |^] ~ [^] . 
Since these motivate the subsequent proposals, the section begins with a 
review of this material. 



2.1 Known cases 

The initial goal is to find sets of equations consistent with movement from 
c = 2(A'^— l)/(A^+2) to c = 1. For even such behaviour has been seen 
before. Each Wso{m) series of coset conformal theories allows for interpo- 
lating flows, analogous to the (/)i3-induced flows in the minimal series. The 
final step, 

50(m)(2) X SO(m)(i) SO{m)^^^ x SO{m)^^^ 

involves a change of central charge 2(2m— l)/(2m+2) 1, which fits the bill 
for N=2m. This was noted by Fateev [^], and our only contribution here 
will be to exploit the TBA systems that have since been proposed: dn+i- 



related for m=2n+2 and 6n,-related for m=2n+l [19, |20|] (n has been 
picked in this way for later convenience: in particular, the choice means that 
the associated sine-Gordon models always have n—1 breathers). The lack of 
an explicit symmetry in the models ( |2.1| ) is not necessarily a problem, 
particularly as attention is being restricted to the flow of the ground-state 
energy. The real test will come later, when a more detailed look is taken at 
the predicted asymptotic behaviours of E{R) at small and large values of R. 
Thus, for A^=4n+4 we conjecture that the ground-state energy will be 



as in ref. |18|: 



re{r) = --Y: / d9ul''\e)4\e) (2.2) 

a = l,2 

where, both here and subsequently. 

The functions el'^\9), i = 1 . . . n+1, are to be found as the solutions to the 
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following set of coupled integral equations: 



n+l 



(a = 1,2), (2.3) 



with Q = 3— Q, and * denoting the convolution 

1 r°° 

f*9i0) = ir / f^^'^9{o - o')de' . 

Zn J-oo 

The kernels and Tpij can be extracted from equation ( |A.8| ) of appendix A 
on setting h=2n. They can equally be defined using the dn+i Toda S- 
matrices, though this fact will only briefly be relevant below. The 'energy 
terms' v\°^\9) are 

z.f ) {9) = {-0) = ^MiRe'^ {i = l... n+l) , (2.4) 



and the numbers Mj are given by ( |A.3| ),( [A^ ), again with h=2n. The full 
system has a ^2 symmetry under a ^ a, 9 ^ —9. This transformation has 
a simple interpretation, as follows. The pseudoenergies e^^' and ef^ can be 
associated with right and left moving massless 'particles' [^], with S-matrix 
elements 

Sf-'^i9) = S^j^i9) = S.,{9) , S^^{9) = S(.^{9) = m,i9))-' 

(with Sij and Tij as in appendix A). The ^2 transformation reverses their 
spatial momenta (by sending 9 to —9) and swaps right and left movers; it 
therefore implements parity. In the rest of the paper, this sort of language 
will often be used, even though the general problem of associating S-matrices 
to the flows is being left for future work. 

Were the theory conformal, E{R) would be given in terms of the central 
charge c by the relation E(R) = —ttc/GR |2^. For a non-conformal theory 
an 'effective central charge' c{R) can still be defined, as 

E{R) = -^c{R), 

after which the ultraviolet and infrared central charges can be extracted as 
the R^O and R^oo limits of c(i?).0 These limits fohow from (U) and (U) 



^Strictly speaking, any bulk terms must also be subtracted before the infrared result 
is valid. The TBA does this automatically. 
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via by-now standard manipulations, and are indeed as claimed. For the 
moment the only point to note is the way that the i?— >oo limit works: 



(2.') 1 

L\ (6) becomes vanishingly small at any values of 6 for which L\'''{9) is 



(1)/ 



non-zero, and vice versa, the two halves of ( |2.3p decouple, and the equations 
become 



n+l 



[a 



1,2). 



(2.5) 



Noting that the /^-dependence can be removed by shifts in ^ of it log R, these 
'kink systems' can be recognised as two copies of the scale-independent TBA 
equations associated with the minimal dn+i-related S-matrices, and also, 
more suggestively, with the sine-Gordon model at the reflectionless points 

Svr _ 327r 

n+l ^ N 



(2.6) 



For A^=4n+2 the story is much the same. The underlying massless 
scattering theory is not expected to be diagonal, and as a consequence the 
TBA system includes some auxiliary functions, in addition to those 
which appear explicitly in the formula for the ground-state energy. These 
functions, often referred to as 'magnonic' pseudoenergies, make themselves 
felt via their appearance in the coupled set of integral equations that all 
the pseudoenergies must satisfy together. For the (6„-related) case in hand, 
three such functions must be introduced. It will be convenient to label them 
as e^' , e^' and e^' , with the remaining non-magnonic pseudoenergies being 
e'f^ and ef\ i = 1 . . . n. (This notation is in line with the idea that it is 
particle number n that should be blamed for the magnonic terms in the 
equations, n being the label of the soliton in the related sine-Gordon model, 
while 1 . . . n— 1 label the diagonally-scattering sine-Gordon breathers.) The 
appropriate 6„-related TBA system can then be written as 



(2.7) 



with 



ct>,,*Lf\e)-i,,,*Lf\e) 



-knjztfct^i^Li^Ke), (a = 1,3); 
/3=0 
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4"^W = -E41Vi*4''H^), (a = 0,2,4). (2. 

(3=0 



Setting 5 = 4— a, the energy terms are again given by equation ( |2.4j ). The 
masses Mj, and the kernels (pij, ipij and </>i, can be found in appendix A, 
with h = 2n—l = N/2—2. Parity symmetry is implemented as a — > 5, 

e -e. 

In the R^oo limit, the interesting behaviour is confined to kink systems 
near 9 = it log i?, and again the equations separate, into one set in which 
the index a takes the values 0, 1, 2 and a second in which it takes the values 
2, 3, 4. As before, the resulting equations could have been found by exam- 
ining the ultraviolet limit of a sine-Gordon TBA system, this time at the 
coupling 

2 _ Stt ^ 327r 
^ " n + 1/2 - N ■ 

Having observed this connection, we can hope to use it in reverse to 
construct new TBA systems for the remaining (N odd) cases, using as an 
input the known sine-Gordon TBA systems at yet further values of the 
coupling f3. The question is, just which values? Looking at equations ( |2.6| ) 
and (|2.g|), a natural guess comes to mind, and to confirm it we now recall 
the one other case where an exact conjecture for the ground-state energy 
has previously appeared, namely the A'^=5 flow from c=8/7 to c=l. This 
was found by Ravanini et al in the course of an exhaustive study of TBA 
systems of a particular type, and reads |21]: 



1 V 1 

RE{R) = — / do v^^\e)L^^\e) + v^^\e)L^^\e) , (2.10) 

27r J -00 L J 

where the functions £'^'^\6) and £'^^\6) appearing in L(i) and L(6) are cou- 
pled, together with four auxiliary functions e^'^^ . . . e^^-* , in the following sys- 
tem of equations: 

eW(0) = z.W(0)-^4^«V*i^^^W (a = 1...6) . (2.11) 

/3=1 

Here z/(°)(0) = ^MRe^i^ + iMRe-^Jga , ^{0) = l/cosh0, and l^^f is the 
incidence matrix of the eg Dynkin diagram, labelled so that 1 and 6 are 
the two extremal nodes. The ^2 symmetry of this diagram, combined with 
6^—9, exchanges left and right movers and thus implements parity. The 
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R ^ oo limit removes all trace of e^^^ from the equations satisfied by e^^\ 
and vice versa; thus in this limit 



5 



= -^^[JV*^^^)^ (a = 1...5), (2.12) 

/3=1 



with l^^p the incidence matrix of the Dynkin diagram. (An analogous 
equation determines the infrared form of e''''(^).) Just as in the earlier 
cases, the final TBA system could equally have emerged in a discussion of 
the thermodynamics of the sine-Gordon model. This time, the coupling to 
take is = 327r/5 |23], precisely the value predicted by equations ( |2.6| ) 
and (I2I 



2.2 New massless TBA systems 



By now it is natural to suppose that a system of equations appropriate 
to the A^=7 flow will have something to do with the sine-Gordon TBA at 
= 327r/7. In its ultraviolet limit, this latter system becomes 



e«(e) 
e(2)(^) 



(^) _ ,^3*(l(4) (0) + {6)) - ,/.4*L(2) (9) 
(/)2*(ir(^)(^) -LW(0)) 



(4) 



where 



L^'") = log(l + e" 



= log(l + e' 



/I) 



(0) 



(2.13) 



and the kernels are given by equation ( A. 12 ) with h = 3/2. (One way to 
obtain this system is to Fourier transform the Y-system given in ref. j24U , 
divide through by appropriate hyperbolic cosines, and then transform back.) 

The scale R enters into these equations in a trivial way, and can be 
removed by a shift in 9. This is as it should be, since the ultraviolet limit 
has already been taken. The scale-dependent sine-Gordon TBA could be 
recovered on redefining vi^\6) = MRcosh9, a manoeuvre that leaves the 
R^O limit unchanged but forces cs'g(c«)=0, as appropriate for a massive 
theory. However, here the idea is different - we want ( 2.13 ) to describe 
the infrared destination of some massless flow, rather than the ultraviolet 
limit of a massive one. To this end, the scale-invariant system, thought of 
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as describing only leftmoving particles, should be augmented with an extra 
piece to describe the right-movers. This new piece will re-introduce a scale 
dependence, but in such a way that the system ( 2.13 ) is recovered as i? ^ oo. 
Furthermore, the augmented system must treat left and right movers on 
an equal footing, with a 7L2 symmetry exchanging them. Remarkably, the 



system ( 2.13 ) is such that both of these requirements can be satisfied in 
a non-trivial way. On noticing the TL^ symmetry already exhibited by the 



mutual couplings of the magnonic pseudoenergies e 



(2) 



it is natural 



to extend this to the full system with the addition of a single right-moving 
pseudoenergy e^^^(Qy. 



^(5) 



r{6) 



(0) 
(0) 



1,(6) _ <^3*(L W {9) + {9)) - (l)4*L^^^ {9) - (l)5*L^^^ (0) 

(2.14) 



^MRe ^ gives the desired left-right symmetry, under e^^^<^e 



Taking i/^^^ 

e(^)<-^e(^), 6^—9. If the ground state energy is 



(6) 



vr 1 /■°° 

RE{R) = --c{R) = -— d9 
b ZiT J -00 



u'^'^\9)L^'^\9) + u'^^\9)L^^\9) 



(2.15) 

then the result c(oo)=l follows, as intended, from the result csg(0)=1 for 
the sine-Gordon model - the determining equations reduce to ( 2.13]) in both 
cases. More interesting is the i?^0 limit of the extended system ( |2.14 ), 
and this confirms that the discussion thus far has been on the right track. 
Using the dilogarithm sum rules given in appendix C, it can be checked that 
c(0) = 4/3, exactly as required. 

Once this case has been understood, the general procedure is clear. 
Take the sine-Gordon TBA at /3^=327r/A^, and look for a ^2 symmetry 
of its magnonic part, which is not a symmetry of the full system. Now 
add pseudoenergies so that this becomes a symmetry of the whole, and 
then check that this enlarged system does indeed give the desired value, 
2(A^— 1) / (A^-|-2), for the ultraviolet central charge. For even the magnonic 
structure is rather trivial and the results are already contained in equa- 

However the idea also works for every odd 



tions (13), dU) and 



N; the resulting systems will now be given. 
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When N=An+l, the magnonic structure of the sine-Gordon TBA is 
identical to that at N=5, as given in equation ( 2.12 ) - the only differences 
come from the appearance of n—1 breathers in the spectrum in addition 
to the fundamental soliton-antisoliton doublet. This leads to the following 
proposal for an A^=4n-|-1 TBA system: 

re{r) = --Y: / deu''\e)L'\e) (2.16) 

27r ^ 7_«3 



with 



(0) 



6 



(Pi,*Lf{e)-i;i,*L'f\e) 



knY. 41^2*4^^^) , («=1,6, a = 7-a) ; 

(3=1 



13=1 



(a = 2... 5). 



(2.17) 



The energy terms are v\ {6) = vl { 



^MiRe^, with the Mi given 



by eqs. (|A.3|) , (|A.4D with h=N/2—2, and the kernels are given by equa- 
tions (ra- dAl^ ). 

Finally, for N=4:n-\-3 the sine-Gordon magnonic structure mimics that 
for N=7, /32=327r/7. The system can therefore be found simply by tacking 
the appropriate number of breather- like pseudoenergies onto equations ( p. 14 ), 
(|1|). With energy terms and kernels extracted from appendix A as before, 
it reads 

RE{R)= / d9 4''\9)L^r\0) (2.18) 

Q = l,6 

with 



- S,,n [MLI^HO) + (0)) + 04*4'^ (^) + 05*4'^ (^) 
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^2^{K!i\e)-L^^){e)) 



n 



(2.19) 



2.3 Modifications for massive flows 

In their ultraviolet limits, MR ^ 1, all of the TBA equations again sim- 
plify, though in a rather different manner to the infrared situation already 
described. In a central region — log(l/Mi?) <C ^ ^ log(l/MR) the energy 
terms v'f^ (6) can be ignored and the pseudocnergies are approximately con- 
stant. Then, in the neighbourhoods of ^=± log(l/MR) , the equations reduce 
to kink systems of a new type, different from those encountered in the in- 
frared in that they continue to involve all of the pseudoenergies. They are 
obtained from the original TBA equations by neglecting either the energy 
terms proportional to MRe^ (this near 0=— log(l/Mi?)), or else the terms 
proportional to MRe~^ (this near 9=+ log(l/Mi?)). If the solutions to these 
two kink systems are denoted and ef respectively, then as i? — 



Once these two kink systems have decoupled, all scale dependence is lost 
from the equations determining the effective central charge, and with the 
aid of the dilogarithm identities listed in appendix C, the value of c(0) can 
be extracted. The result is 2(A^— l)/(A^+2); and this, together with the 
result c((X))=l already found, is exactly as would be expected of a flow from 
the Fateev-Zamolodchikov multicritical point into the Kosterlitz-Thouless 
phase. Subsequent sections will give further evidence for this interpreta- 
tion. But first, the picture should be completed with a set of ansatzc for 
the opposite perturbation, into the massive phase. There is a standard pro- 
cedure to go from a massless to a massive TBA system: simply remove all 
energy terms u^^^ proportional to \MRe~^ , and in the remaining terms vf' 
replace ^MRe^ with Mi? cosh throughout. Up to a relabelling of the nodes 
a— >a for the kink system near ^=— log(l/MR), the leading ultraviolet be- 
haviour of the equations is unchanged by this manoeuvre, and so the result 
c(0) = 2{N—l)/{N+2) is preserved. However, the decoupling previously 




(2.20) 
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induced by the energy terms in the large R hmit no longer occurs, and in- 
stead one finds c(oo) = 0, as required for a massive flow. These systems 
are thus good candidates to describe the situation for the opposite sign of 
the coupling constant, and their properties will also be investigated below. 
It is interesting that after this final step, the sine-Gordon scaff'olding, by 
which we levered ourselves up from c=l to c=2(iV— l)/(A^-t-2), has been 
completely removed. 

Without a prior knowledge of the conformal field theories, the notion to 
reinterpret suitable sine-Gordon TBA systems as the infrared limiting forms 
of larger sets of equations might conceivably have led to their rediscovery. 
It is therefore worth investigating whether there are any other values of 
at which this might work. Curiously, it seems that the points (3^ = 327r/iV 
already captured really are special in this regard - for no other values of the 
coupling constant does the magnonic part of the sine-Gordon TBA system 
exhibit a ^2 symmetry with the desired properties. The only other option 
is to aim for an enlarged system in which the left-right ^2 symmetry acts 
trivially on the magnonic pseudoenergies. Formally, this can be done at 
any value of in the repulsive regime {^"^ > Air; no breathers) it leads 
to sausage models in every case. In the attractive regime the picture 
is not so straightforward. For example, the most natural choice of kernels 
to link the doubled- up nodes, namely the functions ipij successfully used 
above, does not seem to yield any sensible (that is, rational) values for 
the ultraviolet central charges. This is probably related to the 'intolerable' 
properties that the massless sausage S-matrices acquire as the attractive 



regime is entered |25] 



3 Ultraviolet asymptotics 

If the functions E{R) defined in the last section really are the ground-state 
energies of certain perturbed conformal field theories, then their behaviours 
at small R should be consistent with the results of conformal perturbation 
theory. For a unitary theory perturbed by a primary operator e of dimen- 
sions (A, A), this predicts the expansion 

E^^-'\\R) = -^ + '-^Y.Bn.t- . (3.1) 

m=2 

The coupling A appears in the dimensionless parameter t = — 27rA (i?/27r)^~^^ 
while the conveniently-normalised coefficients are given in terms of the 
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connected correlation functions of the perturbing field e on the plane: 



Br, 



1 



(27r)™-im! 



/ m— 1 
{€{l,l)e{zi,zi) . . .e{zm-i,Zm-i))c n 



J, {zuzuY-^ ■ 
(3.2) 

The conformal mapping z = e 2™/^?^ from the tu-cylinder of circumference 
R to the infinite z-plane, is responsible for the appearance of the conformal 



anomaly c in (3.1) and also for the measure used in (p.2|); see for example 
ref. |l^. If 2A > 1, then a finite number of terms, up to order n < 1/(1— A), 
suffer from ultraviolet divergences which must be regularised: the usual 
prescription is to analytically continue in A from a region where the relevant 
integrals converge. Sometimes even this fails to render the answer finite, a 
sign of an irregular term in the expansion of £'(P''''*) . 

The constraints of conformal invariance are such that B2 and are 
completely fixed by the values of A and the operator product coefficient 
Ceee, irrespective of any other details of the theory. For later reference, they 



are |6| 



B, 



1 7(A)2 



a« 7(A/2)= 



4 7(2A) ' ^ 48 7(3A/2) ' ^ ' 

where 7(2;) = r(x)/r(l— x). For the particular models under discussion, the 
values of C^ee were found by Zamolodchikov and Fateev to be S: 



2 n('^'(7)n('^'(2) 



35 



n('v)(4) ^ 



53 



r(#±i)r(i 



y(m±\yi-^ , 5 ^ 



(3.4) 



llfc=i r(i-^) 



where n('-)(j)=nLi ^,._^y 

The perturbative expansion is expected to have a non-zero radius of 
convergence, and the function that it defines should behave at large R as 

being the bulk contribution to the ground-state 
energy. Whilst this is subtracted off in the TBA, this is not the case in 
conformal perturbation theory and £{X) may be nonzero. To take this into 
account, the predictions of the last section should be compared not with 
£;{pert) (lij-gctiy^ rather with the subtracted quantity 



E{X, R) = E^P""''^ (A, R) - £iX)R . 



(3.5) 



Finally, note that the functions E{R) of the last section all in fact depend 
on R and M, where M is a mass scale, part of the infrared specification 
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of the model. (In the massive direction, M will turn out to be the mass 
of a certain asymptotic one-particle state; in the massless direction it sets 
the crossover scale.) Conformal perturbation theory sees instead R and the 
coupling constant A. However, on dimensional grounds there must be a 
relation between A and M of the form 

A = , (3.6) 

where y = 2— 2A is the dimension of A, and k is some dimensionless constant 
relating the short and long distance descriptions of the theory. Later we 
will conjecture an exact expression for k, but for the moment it must be 
considered as another unknown. 

Returning now to the equations of the last section, the first piece of 
analysis, following Zamolodchikov [^], determines the value of y. This relies 
on the fact that any solution to one of the TBA systems above automatically 
furnishes a solution 

to a set of functional relations known as a Y-system. The Y-systems ap- 
propriate to the problem in hand are given in appendix B, and numerically 
they turn out to entail the following periodicity property for their solutions: 

5^/"^ [O + 2vrz^) = y/") (6) . (3.7) 

(a) N-4, 

As a result, the 1^ have an expansion in powers of u{6) = . For 

the full TBA system this will be in positive and negative powers, since in 
the full system energy terms blow up both at 0=—oo (n=0), and at 9=+oo 
(n=oo). However for the kink systems which govern the ultraviolet limit, 
energy terms diverge in one direction only and so = exp e^["^ is regular 

at u=0, and Yi^(") is regular at u=oo. Therefore ef has an expansion in 
non- negative powers of u about u=0, and sf^"'^ an expansion in non-positive 
powers about u=oo. This information is important, because a part of the 
corrections to the asymptotic form of each kink solution can be traced to the 
presence of the other kink at a finite distance in 9, namely 2 log (1 /MR). If 
this is incorporated in an iterative fashion, then the corrections will enter as 
powers of w(2 log(l/MR)). In turn, this implies that in the TBA corrections 
to RE{R) form a regular series: 

oo 

RE{R) = -— + 2^yFm {MRfM^^ , (3.8) 
6 1 

m=l 
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modulo further terms, to be mentioned shortly, that can be traced to the 



subtraction of the bulk term performed in equation (3^). This agrees 
with the result ( p.lD from conformal perturbation theory so long as y = 
2(A^— 4)/(A^+2). This gives a first check on the TBA: the prediction 

6 



N + 2 



for the conformal dimension of the perturbing field. From the work of 
Zamolodchikov and Fateev Q, there is indeed a suitable scalar operator 
in the 2 jy-symmetric conformal theory: sometimes denoted e^'^^ it is both 
Z^v-neutral and self-dual - precisely the two properties that were demanded 
in the later part of section 1. Furthermore, Fateev |17| has shown that this 



perturbation is integrable, so the existence of a TBA system describing it is 
to be expected. 

Next we turn to the bulk piece, which, via equation ( |3.5D , should show 
itself in the ultraviolet as a term £{X)E? in the expansion of RE{R). Fol- 



lowing refs. 1 13, |2^, it is possible to extract the exact value of the coefficient 
£, as a function of Af, from the TBA. The result, which depends on whether 
the system is massive or massless, is 



--massive — 



_ sin6'Arsin26'Ar _ M'^sm^ON 

^massless — 'Z : ^T/i ' W""^/ 



sinS^AT ' 2 sinS^AT 



where Ojsi = 2tt/{N—4:). This formula can be derived generally for A^>8, 
and was also checked by specific calculations for A^=5 and A^=6. When 
A^=4n-|-4, the value of <?massive agrees with an earlier result |17] for the mas- 
sive (/'i,i,Adj perturbation of the d^^^ixd^]_i/d^^^i coset model; this follows 



the correspondence (|2.lD already mentioned for some of the massless flows. 

There are two values of N for which these formulae break down: if A^=7 
or 10, the results are formally infinite. This infinity should cancel against a 
term in the regular expansion (|3.1|), to render the overall result finite |29]. 



For the -R-dependencies to match up, m, the order of the putative cancelling 
term, must be equal to (A^+2)/(A^— 4); this is indeed an integer for A^=7 
and 10. (Such a resonance is also found in the A^=5 theory, between the 
bulk term and the term of order in the regular expansion, but (^]^) shows 
that the divergence is anyway absent in this case.) The dominant residual 
piece can be extracted by evaluating (in the massive case) 



lim 

. 2 2m + l 

m— 1 



2sm39N V ' ^ ' 
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and then setting m = (A^+2)/(A^— 4) at the end. This gives 



, sin 6'Ar sin 26'Ar , 

i^massive = "M^ ^ \og MR . 3.10 

The calculation for the massless case is almost identical: just replace sin2^7v 
by sin^AT throughout. The two instances of these formulae pertinent to the 
current discussion give us 



(7V=7) ^ log Mi? f(^=.10) = AM2 

'massive 47r ^"-"B ) '-'massive 87]- 



i'lsiel = -^M^logMi? , = ^MHogMR 

(3.11) 

It is also possible to say something about the nonsingular parts of the bulk 



terms, which have so far been disregarded: see equation ( 3.14 ) below. 

These results have been obtained by analytic continuation in N , back 
to the 'critical' values of 7 and 10. But it should also be possible to see 
the logarithms directly in conformal perturbation theory, at order when 
N=7 and when A^=10. Note how this is already reflected in the relative 



signs of the massive and massless results in ( 3.11| ). To say more, the inte- 



grals ( p. 21 ) must be re-examined, since their values are needed at precisely 
the points where the general formulae ( |3.3| ) break down. The relevant inte- 
grals acquire logarithmic divergences at short distances, though not in fact 
at long distances: the 'metric factors' l/(zfcZfc)^~^ prevent this. To cope 
with the short-distance divergences, an explicit cutoff a can be introduced, 
with the understanding that when the perturbative integrals are done on 
the cylinder, no two operators should get any closer than a. Transforming 
this onto the plane, via z = ^ to leading order in a operators are now 

required to stay at least 2Tia/R apart. This can be implemented by insert- 
ing Heaviside step functions 9\x\ = ^(1^1 ~^^) ^ appropriate. Since the 
interest is in the limit where all the z^^l, the metric factors can be ignored 
and replaced with an explicit cutoff in the infrared. Then the logarithmic 
divergence in i?2, which occurs when A^=10 and A=l/2, is found by 

B2 = J^'^^\e{0,0)e{z,z))ce\,\d^z + n.s. 



1 



2 J2-Ka/R 



r ^dr + n.s. 
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Here, 'n.s.' denotes the non-singular terms which do not contribute to the 
logarithm. Similarly, when N=7, A=2/3 and there is a logarithmic diver- 
gence in B3. This can be extracted as follows: 

i / |y|<i. 



B3 = J {e{^,O)e{x,x)e{y,y))c0\x\G\y\O\x-y\d xd y + n.s. 



247r2 J [xxyy{x—y){x—y)Yl'^ 



+ n.s. 



' 247r2 J \xx\ 



[yy(l-y)(l-y)]2/3 



d^x + n.s. 



Now the integral in square brackets, call it I{x,a), is convergent (and non- 
zero) at x=a=0, and so to leading order it can be replaced by 1(0, 0), which 
is equal to ttj{1)^ [|^]. The remaining integral gives the logarithm, just as 
in equation ( |3.12D . Collecting the pieces together, 

B3 = C<fr^7 (^)' log ^ + n.s. , (3.13) 



where C^^"''' is given by equation (| 

Since these calculations have left the nonsingular parts of B2 and 
completely uncontrolled, it would appear that no prediction can be made 
about the contribution to RE{R) that is proportional to R^ itself. However 
this is not quite true: assuming that the same regularisation scheme is used 
in the massless and massive directions, all of the perturbative contributions 
*o ^m^siess ^'^'^ ^massive ™ust Cancel out when their sum is taken for N=7, 
or their difference for A^=10. The residual finite pieces then come entirely 
from (|3.9D , and are: 

\/3 

^^massive i "^massless ~ i g ^ ■ (^-l^) 

The plus signs apply when N=7; the minus signs when A^=10. 

Returning to the singular terms, substituting equations ( p. 121 ) and ( p. 13 ) 



into the general expansion ( |3.1| ) gives new expressions for the logarithmically- 
divergent pieces of the bulk terms. But whereas the earlier expressions. 
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equations ( 3.11| ), were in terms of M, this time the couphng constant A 
appears. Equating the alternatives allows the constant k appearing in the 
relation ( |3.6| ) to be determined exactly, at least for these two values of N. 
The results turn out not to depend on whether the perturbations are in the 
massive or the massless directions, and are: 



N=7 : 



N=10 : 



4vr37(i)^G 

3 
8^ 



3 /-y(JV=7) ' 



(3.15) 



There are other values of for which k is known exactly: a special case 
of a formula due to Fateev gives its value for the massive flow, whenever 
N=An ||3^. It turns out that when continued to the fractional values n=7 /A 
and n=5/2, the results ( |3.15| ) are reproduced. This motivates the conjecture 
that the analytic continuation of Fateev's result holds for all values of A'^, 
and for both the massive and the massless perturbations: 



Any N : 



4 7(iv^) 7( JV+2 / 
9^ 



7 



N+2 ) 



( 



r 



2 

JV-4 



JV 
JV-4 



JV+2 



(3.16) 



It is not meant to be obvious that the formulae ( 3.15| ) are special cases of 
this! Further support for the conjecture comes from numerical results, the 
subject of the remainder of this section. 

One reason for undertaking a numerical study of the TBA equations is 
the lack of any analytic results for the coefficients Fm in the expansions 
( |3.8D . In particular, to be sure that the massive and massless TBA systems 
really do describe perturbed conformal theories differing only in the sign 
of a coupling constant A, the relative signs of the F^'s predicted for each 
pair of flows should be examined. Letting Fm denote the coefficients for the 
massless flow, and F^a those for the massive flow, and accepting provisionally 
the conjecture that the mass scales are equal in the massive and massless 
directions, the expectation is 



{-irF„ 



(3.17) 



An iterative method was used to solve the TBA equations for N=5, 
6, 7, 9 and 10, discretising the 9-axis in steps of 69=0.1 and normalising 
1. Direct iteration fails to converge |31]; adding ej"'* to both sides of 



M 
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each equation, dividing by two and basing the iteration on the result solves 
this problem, though convergence is rather slow (up to 300 iterations were 
used at each value of R). Results for RE{R) to about 14 digits precision 
were obtained at 30 values of {MR)y between 0.03 and 0.5. The exact 
constant piece, namely — 7r(A'"— l)/3(A'^+2), was then subtracted, along with 
the bulk terms as predicted by equations (|3.9D and ( |3.1l|) . In addition, for 
N=7 and 10, the unambiguously nonperturbative contributions to the bulk 
terms, identified using equation (3.14), were subtracted in a symmetrical 
way. The remainder of RE{R) after all this should in each case have a 
regular expansion in {MR)^, with coefficients that conform to the general 
rule (|3.17| ). A fit against the numerical data was consistent with this, and 
the values of the first few Fm and Fm thereby obtained are reported in 
appendix D. Note that the TBA system for A^=8 is simply a doubling-up 
of Zamoldchikov's proposal for the tricritical Ising to Ising flow [|^]. The 
expansion coefficients for A^=8 are therefore easily obtained from the earlier 
data, and are also included. In our fits, the m=0 and m=l coefficients were 
left free; comparing their measured values with the exact predictions that 
they be zero gives a check on the numerical accuracy. 

Besides the clear agreement of the results with equation ( p. 17 ), a number 
of other features deserve mention. 

First, to within the numerical errors the ^5 TBA has F^ = F3 = 0. This 
matches with the perturbative expansion: when N=5, 11*^' (7) = and, 
from equation (3.4), C^H vanishes. 





Bl/Bi 


Fi/Fi 


Fi/Fi 


6 


0.046260423 


0.046260427 


0.046260424 


8 


0.379827746 


0.379827746 


0.379827746 


9 


14.448969 


14.448967 


14.448991 



Table 1: Comparison of the Zg, and systems with conformal pertur- 
bation theory 

For A^=6, 8 and 9 a more delicate test can be done. Independently of 
the value of k, F^/F^ = F^/F^ = Bl/Bj should hold if the TBA is to 
be consistent with conformal perturbation theory. Using equations (|3.3| ) 
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and ( |3.4| ) to evaluate iJf/^l, the numerical results are indeed consistent 
with this, as table || demonstrates. 

Finally, the results for N=7 and 10 support the formulae (3.11) and 
( |3.14 ): any errors in the bulk terms subtracted would have destroyed the 
good agreement that the fits show with the prediction ( 3.17| ). 



N 


Af2 

"'exact 


;f2 

"'num 


num 


5 


0.0235204664 


0.0235204667 


0.0235204663 


6 


0.0371477546 


0.0371477547 


0.0371477547 


7 


0.0434800500 


0.0434800501 


0.0434800501 


8 


0.0442207130 


0.0442207131 


0.0442207131 


9 


0.0418094000 


0.0418094001 


0.0418094003 


10 


0.0379954439 


0.0379954451 


0.0379954487 



Table 2: Comparison of exact and numerical values for 

Once the consistency of the data with conformal perturbation theory has 
been established, numerical estimates for k can be extracted and checked 
against the conjecture ( |3.16D . Equating the terms proportional to B^^ in 
equations (|3.1D and (|3.8D yields 



F2 _ {2T:fy 7 (1^2) 



«num - l^Vrj - ^ t2 , 

and similarly for /5num- For A^=10, the second-order (m=2) term is obscured 
by the logarithm, but a prediction for k can be obtained in an analogous 
way from the third-order behaviour. The numbers obtained in this way are 
compared with the exact predictions, K^^act' ™ table the agreement is 
very good. 
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4 Infrared limits 



The behaviours of the massless and massive systems become very different as 
R grows, and they will be described separately. We start with the massless 
flows. 



4.1 Into the Kosterlitz-Thouless phase 

The massless systems all exhibit c(oo)=l, a prerequisite if the picture out- 
lined at the end of section 1 is to be confirmed. For a more detailed check, 
some information on the possible operator contents of each 2 at model in its 
Kosterlitz-Thouless phase is needed. A complete description can be found in 
ref. 1 32]: the scaling limit of the model is equivalent to a free boson compact- 
ified on a circle of some radius r. Adopting the normalisation traditional in 
conformal field theory, the action can be taken as 

= ^ J dMd.^f . (4.1) 

The boson <I> can be split into holomorphic and antiholomorphic components 
by writing z = xi+iX2 , ^{z,z) = ^{(j){z)+(j){z)) , with 

{(l){z)()){w)) = - log{z-w) , {(l){z)(j)(w)) = - log{z-w) . 

Compactifying onto a circle of radius r amounts to the identification $ = 
$-|-27rr. Closer to the discussion of section 1 would have been to work with 
9{z,z) = (^{z,'z)/r, whereupon the prefactor in ( [4.1D would have become 
r^/27r. This shows that r^/vr should be thought of as the inverse tempera- 
ture. 

Invariance under $ — > ^+2'nr and mutual locality restrict the possible 
vertex operator primary fields to 

ynm{z,z) = V2cos{p<t){z)+p'4){z)) , V~^{z,z) = V2s\n{p(t){z)+p'4>{z)) , 

where (p,|)) = {^+mr,^—mr) and n,m £ '2.. These have the conformal 
weights 

(A„^, A„„) = (i/, ip2) = (i(|.+^r)2, i(|-_rnr)2) . 

Included in this collection are the the vortex fields V^^j, and the A^-fold 
symmetry breaking fields V^q. Their conformal weights are: 

AAro(r-) = AAro(r) = , AoM(r-) = AoM(r) = — — . 
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This information allows the of values of r actually realised in the Kosterlitz- 
Thouless region to be determined. At the low-temperature (large-r) border, 
the fields V^^ become relevant, so Aivo('^max) = 1> and rmax = N/\/8. 
Similarly the high-temperature boundary is encountered when the least ir- 
relevant vortex fields, Fq^, become relevant, that is when Aoi(r) = 1, so 
'"mill = \/2- One further value of r can be identified in the thermodynamic 
phase space. The duality transformation for the rnodel exchanges the 
A^-fold symmetry-breaking field with the unit vortex, that is Vtvo with Vqi. 
Thus Tsd, the compactification radius appropriate to any self-dual system, 
must satisfy Ajvo('^sd) = Aoi(rs(j). Hence, r^d = \/N/2. To summarise: 



r 


Aatq 


Aoi 


''mill — \/2 


jv2 
16 


1 




JV 
4 


JV 
4 


r — JV 
'max— 71 


1 


JV2 
16 



Since initial (ultraviolet) perturbation was self-dual, it is the middle row of 
the table that will be important. The discussion has been couched in the 
language of continuum-valued spins, and it might appear surprising that 
a model with discrete spins can mimic such behaviour. Universality and a 
suitable renormalisation group should anyway suffice, but for a more physical 
picture of the mechanism involved, see ref. p^]. 

Pure c=l behaviour is only expected in the very far infrared. For any 
finite value of R, there will be corrections to scaling stemming from irrelevant 
operators compatible with the symmetry of the model [^]. Apart from TT 
and its descendants, the most relevant of these when the theory is self-dual 
have conformal dimensions {N/4:, N/4). The model is therefore described 
by the continuum action 

S = S'fR + f^ij d'^x + ^^y" TTcfx + (further terms) , (4.2) 



where 5j*j^ is the action (|4.1|) , and 'i/'7V/4 denotes an equal (self-dual) combi- 
nation of the unit vortex and A^-fold symmetry-breaking fields. The further 
terms are not optional in the continuum theory: the perturbation theory is 
non-renormalisable, and so infinitely many counterterms are needed. Never- 
theless, it is possible to extract some non-trivial predictions: an asymptotic 
expansion for RE{R), in various powers of R and /J^i, ■ ■ ■ , can be obtained 



after a mapping from the cylinder to the plane [15 
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The fact that is primary means that its contributions only start 

at order /if , and the charge conservation selection rule for vertex operator 
correlation functions restricts subsequent terms to even powers. For this 
operator, y = (2— 2A) = 2— N/2, and so its leading correction to RE{R) is 
generally proportional to ^\R'^~^ . 

By contrast, TT is a secondary operator and has a non-vanishing one- 
point function on the cylinder. Consequently, despite its renormalisation 
group eigenvalue y being equal to 4, the corrections due to TT start early, 
with a term proportional to ^2R~'^- 

Evidence for 'i/'Af/4 can be found in the influence of the magnonic pseu- 
doenergies on the large R solutions of the massless TBA equations, using a 
small extension of an argument given by Zamolodchikov |jT^. In the R^oo 
limit, the pseudoenergies acquire the form of a pair of sine-Gordon kink sys- 
tems, one located near 9=— log MR, and one near 9=log MR. In the central 
region between the kinks, the non-magnonic pseudoenergies are dominated 
by their energy terms, and play no direct role in the equations: the func- 
tions -L*j^' and L'^\ which would otherwise couple to the reduced system of 
magnonic pseudoenergies, are doubly-exponentially suppressed. Their dom- 
inant effect comes instead via the 'tails' of the relevant kernel functions, 
which only suffer a single-exponential decay. 

Thus for N = An+2 (when the TBA system applies) the magnonic 
pseudoenergy ej^^'(^) sees e^^'(^) and e^^'(^) via the convolutions (^i*L^^'(0) 
and (/>i*L^'(0). In the central region these expand as odd series in e^^^ 

(recall that (l)i{6)=h/ cosh. hO ), implying that Y^\0)=e^'^ satisfies the 
equation 

YP{e-f^)YP{9 + §^) = l. (4.3) 

(Note that this is just the truncation of the full Y-system to the central 
magnonic node.) Therefore ^n^'(^+^) = ^n\^)- Repeating the argu- 
ment used in the ultraviolet limit, this periodicity implies that the leading 
'magnonic' corrections to RE[R) are powers of {MR)~'^^ = {MR)^~^ . 

For = 4n-|-l and A^ = 4n-|-3, the magnonic pseudoenergies form a 
d4-type system with kernel </>2(^) = 2h/cos\i2h0. The pseudoenergies e[i^'(0) 
and s''^\9) are felt via (/)2*L^^'(^) and (/)2*L^®'(0). These expand in the central 
region as odd powers of e^^'*^, and such terms cancel out when shifted by 
ibz7r/4/i and summed. As a result they do not impede the construction of 
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the Y-systems. For = 4n+l, the system is 



yt\o - ^) Yt\e + n) = n (i + Y^ie)) , (4.4) 

while for = 4n+3 , 



13=2 



5 _;[d4l 

yt\0-vi)yt\o + ^^) = \{{i+Yi,^\e)-') . (4.5) 

/3=2 

The index a runs from 2 to 5 in both cases. Again, these are just truncations 
of the fuh Y-systems. They imply the periodicity y^")(6l+^) =Yj^"\d) |2^, 



after which the argument runs just as for A^ = 4n+2, and again predicts 
corrections to scaling as powers of (Mi?)^~^. 

Finally, to the case A^=4n-|-4. There are no magnonic pseudoenergies 
in the TBA system, and so all corrections to scaling must come from the 
'direct' interaction between the non-magnonic pseudoenergies, to be dis- 
cussed shortly. This would seem to be a problem, since the general idea 
will be to attribute these corrections to TT and descendants. Two points 
can be made. First, contributions from ipN/4, expected to appear as powers 
of {MR)~'^^, will anyway be badly tangled up with those from TT and its 
descendants. However this point applies equally when A^ = 4n+2, where at 
least for A^ = 6 the entanglement manifests itself in a logarithm at order 
(MR)'^ (see below). Such a logarithm is not found at order {MR)~^ in the 
^8 theory: in fact the coefficient of this term fits perfectly the hypothesis 
of a pure TT perturbation. (Note, coefficients up to high order can be ob- 
tained for this case by doubling the (exact) results for the tricritical Ising 
to Ising flow obtained by Zamolodchikov |jl6[.) The resolution comes on re- 
examining the formula ( ^.3| ) for the coefficients in conformal perturbation 
theory: when A^=4n, A=n and is a positive integer, and B2 vanishes. Thus 
signs of tpN/i anyway appear anomalously late in the series. Without any 
control of the counterterms, it is not possible to be any more precise than 
this, but at least the proposal has not turned out to be inconsistent. 

So much for For TT, the analysis can be pushed a little further. 

To eliminate the corrections already discussed, it is convenient to replace the 
TBA equations with a set in which the magnonic pseudoenergies associated 
with the two kink systems are separated, by doubling them up and cou- 
pling one copy only to and the other only to e^^K Since the magnonic 
interaction only cuts in at order {MR)^~^ , the approximation should be 
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good at least up to this point. For A^=4n-|-4, no approximation is needed 
at all. The kinks near 0=— log MR and 0=log MR now only interact via 
the tails of the kernel functions tpij{6). In Zamolodchikov's original paper 



on the massless TBA an iterative procedure was used to account for 
such an interaction. It turns out that in all of the subsequently- introduced 
TBA systems for the flows ff^^^ x^^i)/^^^) ^ c/^^) xc/(i)/g(2)^ g e ade 



the first two corrections to infrared scaling also have a simple form. Only 
the case g = dn+i will be directly relevant here, but the general result will 
be recorded as it does not seem to have been given elsewhere. The TBA 
systems were given in their full generality in ref. [^], and involve kernel 
functions 4>ij(P) and ipij{9), related to the corresponding purely elastic S- 
matrix elements in the way described for the dn+i case after equation ( [2.31 ). 
Using some tricks previously employed in the (ultraviolet) calculation of the 
bulk terms for certain massive flows from the g^^^ x g^^^ / g^'^^ cosets [^], the 
final result is 

^REiR) = -^c(oo) + ^ + ^ + 0{{MRr') (4.6) 

where c(oo) is the central charge of the infrared limit. The coefficients Ci 
and C2 are 

C, = -H(^^^41^ , C. = -Hi^f^^4l^V (4.7) 
12 3 ' 6 \^ 3 Mf y ^ ^ 

where Mi is the mass of the lightest particle in the theory, and t/j^^^ is the 
first coefficient of an expansion of the kernel ipij{9) ■ 

00 

s=l 

These coefficients are only non-zero when s taken modulo h is an exponent 
of the relevant non-affine algebra g. An explicit expression is 

/{*) ^ W i^) (A o\ 

h 

where qf^ and q^^^ are components of a unit-normalised eigenvector of the 
Cartan matrix of g with eigenvalue 2— 2 cos ^s. This can be derived from 



the general S-matrix formulae of ref. [36|, using the identity 



h-\ hq^'\^^^ 
5:(A,,i.-^'0,)sin(2p+l+n.,)f = -^—^ (4.9) 

„ Z bill L, 
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which is the Fourier transform of equation (2.14) of ref. ||3^. (See refs. 



37 1 for an explanation of the notation used in this last formula.) 



The results (4.6),( p77| ) can be compared with the perturbative expansion 
of the action S^^ + H2J TT(fx |K 




1 RPrp^ 1 ( ^^(2v^)V2 /c(oo)y (27r)Vi 

' (4.10) 

This much of the expansion is uncontaminated by counterterms, and there- 
fore gives the unambiguous prediction CI/C2 = -c(oo)/24 if the TBA is to 
be compatible with conformal perturbation theory. It is easy to verify that 



this does indeed hold when the coefficients are given by (O). (This was 
checked numerically by Martins, for (7=02, in ref. [^.) A comparison of 
the two series also allows the renormalised coupling (j,2 to be expressed in 
terms of the crossover scale M: 

Note that for any g £ ade, this ^2 is negative: for a unitary flow with TT 
contributing the leading infrared correction, this has to be the case, since 
otherwise the expansion ( [4.1[1| ) would contradict Zamolodchikov's c-theorem. 

Specialising now to g=dn+i, we obtain results relevant to the self-dual 
Zjv flows with A^=4n+4. For other values of A^, the kernel i/jnniG) is more 
complicated. Nevertheless, the identity 

4;' = ^*!;' (U = l...n) 

(previously a simple consequence of equation (|4.8| )) continues to hold and 
allows the calculation to go through as before, so that the results ( |4.6[ ),(4.7) 
are again recovered. Substituting c(oo) = 1, Mi = 2M sin j^^^ , and V'ii = 
4 sin , the final values for the expansion coefficients and TT couplings 
read 



3sin / ' 6 \3sin 



7V-4/ " Y 00111 ^_4/ 

^2 = - , M'^ . (4.12) 

^ sm 
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These formulae break down when A^=5 and A^=6. This is not too sur- 
prising: the whole derivation fails for N<7, since the direct interaction 
kernels tpij{0) are zero in these cases. As the numerical results below show, 
corrections do nonetheless appear. At a qualitative level, these can 
be seen to come, exceptionally, from the magnonic kernels (though for A'^=7 
it is necessary to invoke the thus-far disregarded kernels 4>s, 4>4 and 05 ). 
However, for exact results it is more fruitful to suppose that an analytic 
continuation of the formulae valid for larger values of makes sense. For 
N=7 this is easily done, and gives a result that checks well with the numeri- 
cal data. For A^=5 and A^=6, a little more work is needed: in the same spirit 
as for the bulk term logarithms, the 'TT' divergences given by should 
balance against divergences in other parts of the expansion, allowing the 
final result to remain finite. Indeed, at these two values of A^ an additional 
contribution proportional to comes from the term of order p^^^ in 
the perturbative expansion, and this has the potential to cancel the pole in 
Ci (furthermore, the regularised B2 itself diverges when A^=6). To find the 
residual finite term, set m = 2/[N—4) and evaluate 



lim 



TT / 1 1 



12 3sin^ \(MRY (iVffi)(^-4)'^ 



(A^-4) log MR 



36 cos {MRf ' 
(4.13) 

This gives the corrected coefficients 

C^^=^) = J_ log MR + const , c[^=^^ = _ J_ log MR + const . (4.14) 
36 18 

These values can also be extracted directly from the TEA, as explained in 
section 4.2 below. To deal with C2, two additional terms must be taken into 
account, at perturbative orders ^f"^ and //™/i2 : both become proportional 
to R~^ when A" = 5 or 6. In terms of x = A^ — (4-|-^), they should cancel 
the X— >0 divergences in 



Co(x) 



1 m ?>m? m'^ir'^ 



27m4 

To achieve this, a(x) and b{x) must be found such that 

1 f 1 ^ mx a{x) b{x 



A+mx 



x2 V(Mi?)4 (Mi?)4 (Mfi)4+2mx (JVffl) 

is finite at x=0, determining a(x) = — 1 -|-qx-|-0(x^) , b{x) = 2 + {m—a)x + 
O(x^). In turn this fixes that part of the finite residue proportional to 
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log^ MR. The upshot is that the earher expression for C2 should be modified 
to 

^(N=5,6) ^ iog2 MR + const, log MR + const' . (4.15) 

Infrared logarithms appearing to a single power have previously been ob- 
served numerically, in the interpolating flows between the low-lying minimal 



models Aip Mp-i ||2^, [S^], and also analytically, in a massless sausage 
model p5| . The method used in section 4.2 to extract the logarithms 
directly from the TBA can also be applied to these other flows, and the 
calculations are described in second half of that section. 

We also investigated the infrared behaviour of the TBA systems numeri- 
cally, both to verify the results derived above and to provide information 
on coefficients that we have not been able to extract exactly. For R ranging 
from 120 to 11000, ^RE{R) was found for iV = 5, 6 and 7. The iterative 
procedure was as described in the last section, and again the normalisation 
M=l was adopted. 

First, we estimated the exponents of the leading correction terms by 
finding the limiting slopes of plots of \og{^RE{R) + ^) against log R, with 
the results —0.9993, —1.9967 and —2.0010 respectively (to be compared 
with the predictions of —1, —2 and —2). Thus reassured that at least the 
leading behaviour was as expected, we then fitted the data to expansions in 
R and log R of the predicted forms, leaving all the coefficients unconstrained. 
Numerical precision was not as good as in the ultraviolet, possibly reflecting 
the asymptotic nature of the infrared expansions. The following fits were 
obtained for ^RE{R) + ^ : 

0.0177380 0.027733 log i? 0.01951 0.097 

■ R—^ m ^ + ^ + --- 

^ 0.05555546 log 0.033585 0.0715 log^i? 

■ R- + — R^ + ••• 

^ 0.1007662 0.153 
N=l : ^ + ^ + --- (4.16) 

In each case, the omitted constant piece (measuring the difference between 
the numerical and exact values for c(oo)/12) was zero to 11 significant fig- 
ures. For A^=5 and 6, the coefficients of R~^ log R match well with the pre- 
dicted values of 0.02777. . . and —0.05555. . . respectively, while for A^=7 the 
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formula ( p^ predicts the coefficient of R'"^ to be -0.1007663. . . , again m 
good agreement with the numerical results. The other predictions are more 
difficult to verify: if the expansions really are asymptotic, then results at 
larger and larger values of R are needed before higher terms can be captured, 
and one rapidly runs up against the limitations in numerical accuracy. For 
N=7, a coefficient of 0.24369. . . is predicted for the R~^ term, while for 
N=5 and iV=6 the coefficients of iJ^^log^ R should be -0.0185185. . . and 
—0.074074. . . respectively. The agreement for A^=6 is reasonable, and im- 
proved when a restricted fit (using the exact values for c(oo) and Ci) was 
performed, the estimate for this coefficient changing to —0.0743 .... For the 
other cases, no such convergence was observed for any coefficients beyond 
those reported above. This is not particularly surprising, since for A^=5 
there are more unknowns to fix before the term of interest is reached, and 
for N=7 it is anyway smaller, lacking the log^ R factor. Greater numerical 
precision will be needed before any more can be said. 

Nevertheless, in all instances where we have reliable numerical results, 
they agree with the exact predictions. Since some of these latter were de- 
rived using an additional physical assumption of analyticity in A^, this lends 
further credibility to the whole scenario. 

One further remark can be made before leaving the Kosterlitz-Thouless 
phase. If the action defined by (^]^) and ( [4.1[ ) is converted into the nor- 
malisations of appendix A by setting <I> = y/nip, then the A^-fold symmetry 
breaking part of the perturbation becomes proportional to cos(\/2A^vr (p), 
and formally matches the sine-Gordon Lagrangian (AT) at = 2A^7r. 
Toda-type duality (see for example ref. 09]) maps /32 to 647r2//32 in these 
normalisations. Thus = 2A^7r is Toda-dual to = 32-k/N, the value 
already picked out in section 2. This gives an alternative perspective on 
these special values of the sine-Gordon coupling constant. 



4.2 Exact infrared information from the TBA 

This section, something of an interlude from the main development of the 
paper, discusses some general properties of massless TBA systems. Moti- 
vated by the results obtained above, the particular interest is the mechanism 
by which logarithmic terms can be generated at large R. The simplest ex- 
ample of this phenomenon is the flow between the minimal models A^s and 
A^4 The relevant TBA system, proposed by Zamolodchikov |l^, reads 

iMiJe^ = ei{9) + ^*L2{e) (4.17) 
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= 62(0) + (l)*Li{e) + (l)*L3{e) (4.18) 
^MRe-^ = e^{e) + (p*L2{e) (4.19) 



with (l){6) = l/cosh0, Li = log(l + e ^^), and 

111 1"^ 
—RE{R) = -—c{R) = / deMRe^Li{9) 



oo 
oo 



--^ / deMRe-^L^ie). 

47r^ J_oo 



Take the derivative with respect to 9 of (4.17), multiply by —■^Li{9), and 
integrate from — oo to oo to find 

1 1 roo 1 roo 

—RE{R) = -^j d9 deei{9)Li{9) - ^ J d9 {^*deL2{9)) Li{9) . 

(4.20) 

At R=oo, El solves a kink system obtained from equations ( 4.17D -( |4.19D 



by omitting all terms involving £3, and vice versa for £3. In line with the 
notation used earlier in the paper, write ef , solutions to the 

first kink system, and efi ^3 ^r the solutions to the second. To find the 
logarithmic correction, which comes from an integral over the central region 
— log MR <C ^' <C log MR between the two sets of kinks, replace Li with Lf 
and d0L2 by SgLf + dgL^- This gives 

1 1 1 

l^RE{R) = - Y^c(oo) - ^ jj9 {<^*deLm) Li{0) ■ (4.21) 
In the central region e^^') satisfies 



1 1'°° 1 

J_ / ^n/f j_ 

2tt 7_oo cosh(e' - 9") 



= £^{9') + —l d9" — — — —Ll{ 



p" /-oo 

eU9') + — d9"e-'LU0") + 



IT 



00 



and so 



- (4.2.) 

,4.23, 

Because of the double-exponential suppression of Lf{9) in the central re- 
gion, (j}{9-9')deL2{9') can be replaced there by 2e^~^' deL2{9'). It therefore 
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has a piece constant in 0', equal to e^7rc(oo)/3-/Vffi. Dividing by 27r and inte- 
grating across the central region for the convolution then produces a factor 
-log MR, and equation (4.21) becomes 



—RE{R) = -— c{oo)- ^ r dOe^^^ log MR Lf (9) + .. . 
2iT ^ ^ 12 ^ ^ 27r2 i„oo 3MR ^ ^ ^ 

1 c(oo)'^ 

With c(oo) = 7/10, the coefficient of the logarithmic term is —49/1800 = 
—0.027222..., which compares well with the value —0.02723(2) found nu- 
merically by Klassen and Melzer in ref. [p^ j. 

The calculation goes through essentially unchanged if the TBA system 
is augmented by equal numbers of rightmoving and leftmoving pseudoener- 
gies, so long as they interact only with ei and £3 respectively. Their effect 



is to change the value of c(oo), and the result ( 4.24 ) remains valid. In par- 
ticular, the W.Q TBA system, a special case of ( p.8[ ), can be obtained in this 
way, simply by adding two magnonic pseudoenergies eo and £4 to the sys- 
tem ( [4.17 -4.19). (Note, when A^=6, (j)i{0) is exceptionally equal to l/cosh0.) 
Since c(oo)=l for this system, the value of C}'^^^' given by equation ( 4.14 ) 



follows immediately. 

A different generalisation is needed to find Cj^^""', since for A^=5 the left 
and right kink systems influence each other via a reduced system containing 
more than one magnonic pseudoenergy. This will modify the result ( 4. 231) , 



but if dgL2{9') retains a piece proportional to e^' in the central region, a 
logarithm will still appear. Formally the question is the same as arises in 
Zamolodchikov's treatment of ultraviolet logarithmic terms, described in 
appendix A of ref. |1G], since to leading order the equations are the same 



as those for a standard TBA. The one change is that instead of energy 
terms ^MRe^^, the reduced system is driven by terms ^^mr e^^, appear- 
ing for those pseudoenergies which in the full system would be coupled to 
the remaining right or left moving pseudoenergies. These are induced in 
the same manner as already seen in equation ( [4.221 ). (An analogous effect 
was important in the derivation of equations ([4.3[ )~(^^) of section 4.) By 
reversing some of the steps of Zamolodchikov's ultraviolet argument, and 
then repeating them in the new context, it is possible to derive the following 
prescription. Suppose that the TBA based on the reduced set of pseudoen- 
ergies would have shown an ultraviolet logarithmic correction to ^RE{R) 
equal to F'^^^ {MR /I-k)"^ log MR, had the energy terms been of the usual 
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form. Then towards the infrared Umit of the full problem, 50^2 (^') develops 
a piece proportional to in the central region: 

dem') = ... + Fl^JJ^e'' + ... {-logMR « 9' « logMi?) . 

Note, in contrast to equation ( [4. 231 ), this will not necessarily be the leading 
term. Nevertheless it gives the leading logarithmic contribution, feeding 
through the remainder of the calculation to predict 

lREiR) = -^^cioo) + ... + ^,lo,MR + ..., 

with Ci = ^i^J^g;c(oo)2 . (4.25) 

The earlier result ( [4.24| ) can be recovered on noting that the reduced TBA 
there is the same as that of the thermally-perturbed Ising model, for which 
^(log) ~ For the ^5 TBA, the reduced system is based on the d4 

Dynkin diagram, with driving terms attached one to each node of the fork, 
and -^(*iog)' = 1/4 |Q . Hence Cj^"^' = 1/36 , as previously obtained by other 
means. 

The same method can be applied to the sausage models SST^ ^ at A = 



1/K, K = 3,4, . . . , using the TBA systems given in ref. |25|. These models 
arrive at c(oo)=l, attracted in the infrared by an operator with conformal 
dimension A//? = 1/(1— A) = K/{K—1). The reduced system is based on 
the (non-affine) dx-i Dynkin diagram, with (for K > 4) both left and right 
moving driving terms coupling to the extremal node of the tail. This entails 
Fl^^^ = for even, and F^^^'^j = -2/{K-l) for K odd Q. Hence 

f {K even) 

When K=2, the reduced system is disconnected, and there are two contribu- 
tions to the logarithmic term. When these are added together, the general 



formula turns out to be valid in this case as well. In ref. |25] the same result 
for K='i was derived by an indirect route, closer to the analytic continuation 
argument used elsewhere in this paper. Part of the input was the pair of 
formulae 



Cf' = -^tan^ , 4^) = -^tan|^M-, (4.27) 
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found by putting the model in an external magnetic field. In fact this infor- 
mation is enough to predict the coefficient of the logarithm for all integer K, 
with agreement with the direct TBA derivation. This gives a further check 
on the results of ref. ||2^. Going further (and contrary to equation (5.47) 
of [p^]), arguments similar to those used in section 4 predict an additional 
log^ffi/(Affi)^ correction for K odd, with coefficient -32/27{K-l)'^ . This 
was checked numerically for K=3. 

One other special case of the general result ( [4.25| ) is of some interest. 



This covers the TBA systems proposed by Zamolodchikov in ref. [41| to 
describe the flows 

MA'^+\k,l) : M{k,l) ^ M{k,l-k) (4.28) 

from the coset models Ai{k,l) = o]^"^ x af^ / a'^^^\ The pseudoenergies live 
on an a^+fc-i Dynkin diagram, with a rightmoving energy term ^MRe^ 
on the k^^ node, and a symmetrically-placed leftmoving term ^MRe~^ on 
the l^^ node. The reduced system operating in the central region sees 
the ai^k-i Dynkin diagram, with the induced energy terms sitting one at 
each end. If l—k—1 is odd then the ultraviolet version of such a system 
shows a logarithmic singularity, with F^^^"^^ = 2(-l)('-'=)/V(^-A;+2) 
Therefore the flow ( [4.28 ) will show a logarithm in the infrared whenever 



l—k is even. The coefficient is given by ( 4.25 ), with c(oo) = c{k,l—k) 



3k{l-k){l+A)/{k+2){l+2){l-k+2) . Setting k=2, 1=4 recovers the re- 
sult, while the cases k=l, l=p—2 pertain to the flows M-Ap^^: M.p —>■ M.p-i 
between minimal models. Therefore, for these flows 

{p even) 

Ci' = { , .^f.-3V2 2c(oo)2 (4.29) 




9(p-l) 



(p odd) 



where c(oo) = 1 — 6/p(p— 1). An analogous calculation for the massless 
flows Hp^2- —>■ ■Mp-i introduced in ref. |Q (where Zp_2 denotes 



the 2p_2-symmetric conformal field theory) finds the same formula ( 4.29]) 



save for the replacement of the prefactor (— l)*-^ ^^''^ by —1. The similarity 
of the result to the sausage formula ( [4.26| ) is not coincidental. As 

recently stressed in ref. |p8|, it appears to be consistent to regard the theory 
Hp_2 as a reduction of the sausage model SSTyJ. Furthermore, in such a 
case the ground-state energy E^{R) after reduction should be equal to that 
of a suitable excited state E^{R) in the unreduced model (this is clearly 
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explained for the reductions of the sine-Gordon model in ref. ^^). For the 
reduced ultraviolet and infrared central charges to emerge correctly from 
this manoeuvre, the excited state must have ultraviolet and infrared scaling 
dimensions 

2Afc(0) = ^ (c'^(O) - c'^m , 2Afc(oo) = ^ (c"(oo) - c'-(oo)) , 

with c"(0), 0*^(00) and c^{0), c^{oo) the ultraviolet and infrared central 
charges before and after reduction. In infrared perturbation theory, this 
is reflected in the modification to the TT corrections to scaling for an ex- 
cited state with infrared scaling dimension 2Afc(oo) [^8|: the result ( [4.10| ) 
still holds, but with c"(oo) replaced by c"(cxd)— 24Afc(oo) = c''(oo) . Com- 
paring the TT expansions for ^RE^{R) and ^RE^{R) , this implies that 
the coefficient of (MR)^'^ for the excited state can be obtained simply by 
multiplying the coefficient for the ground state, Cf, by (c''(oo)/c"(oo))^. So 
long as the mass scale remains the same after reduction, this means that 

A similar result relates C2 to C2 . Resonances may ocassionally cause these 
coefficients to diverge. However, so long as the resulting logarithms can 
be obtained by analytic continuation from non-resonant points, their coeffi- 
cients must be related in the same way: 



2 



C; = ( 31^ I , (4^31) 



This explains the connection observed above between the logarithmic singu- 
larities for SSTy^ and when p is odd. Furthermore, combining ( [4.27| ) 

with ( [4.30| ) yields a prediction for when p is even, namely C^^^ = 0. 

The flow MA'^^ can also be found as a reduction of a theory that flows 
in the infrared to c=l, this time an imaginary coupled sine-Gordon model 
(ISG) 1^^. The flow to take leaves c=l by a relevant operator with con- 
formal dimension Ajjv = wanders around in a non- unitary way, 
and then returns to c=l via an irrelevant operator of conformal dimension 
^IR = p/{p—^) (as usual, there will also be a collection of counterterms) . 
Such models can be defined for all p>2. Unfortunately, we do not know 
very much about their infrared behaviour, and so we resort instead to the 
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following trick, which finally brings the discussion back to the TBA sys- 
tems introduced in section 2. First, parametrise the arriving dimension 
A7v=-/V/4 for the self-dual 2jv flow as Ajv = 1), and then rewrite the 



formula (|4.12 ) in terms of p: Ci^^ = — 7r/(36 sin ^(p— 1)) . This is singular 
when p is an odd integer - precisely the values where the corresponding 
coefficient in the ISG should diverge if its reduction is to reproduce the be- 
haviour of MAp^^ . The coefficient of the logarithm arising at these points 



is, from equation ( 4.13| ), equal to (— 1)^^ ^^/^/9(p— 1) . Comparing with the 



result (|4.29| ) for M-A^\ p odd, and recalling the relation ( 4.31 ), leads to 



the conjecture that the coefficient Ci for the massless sine-Gordon model is 
minus twice the analytic continuation of the 2jv value: 

^(ISG) ^ TT 

^ 18sinf(p-l) ■ ^ ' ' 

We have only checked this directly for one case: in ref. [|^, a TBA system 
was given for the p=2 ISG, which arrives at c=l along operators with the 
same dimensions as those for the self-dual flow. From this system, it is 
possible to extract Ci = vr/lS, which is indeed the value given by ( [4.32] ) . 
Note that the expected monstrous corrections to the ISG TBA, investigated 



in ref. [44|, have an exponentially small influence in the far infrared and so 



do not change this result. At other values of p, the conjecture leads to a 
prediction which can at least be checked numerically: for AiA^^^^ with p 
even, a case previously inaccessible to us, we expect that 

For p=4 this reproduces the value = — it /72 obtained by Zamolodchikov 
for the tricritical Ising to Ising flow [16|. To go further, we extracted values 



of Ci or Ci from fits to numerical solutions of the relevant TBA systems, 
allowing for the possibility of log^ MR/{MR)^ terms for p odd, for p = 
5 ... 10. The results are compared with the exact predictions in tables ^ 
and ^. (The first numerical entry in table ^ was taken from ref. [p8|.) 

The paper |Q, by Feverati et al, appeared as this paper was being 
written. Amongst other things, it contains a collection of further numerical 
results, in particular for the H^_2 flows, which are also consistent with the 
exact results obtained in this section. 
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exact 


numeric 


M5 M4 


-0.027222. . . 


-0.02723 


M^ Mq 


0.027211... 


0.027211 




-0.023341 . . . 


-0.0232 



Table 3: Ci for the flows MAy'^ p odd. 



Mp Mp_i 


exact 


numeric 


Me ^ Ms 


0.111701... 


0.11166 


Ms My 


-0.139137... 


-0.13910 


Mio ^ Mg 


0.152038... 


0.151 



Table 4: Ci for the flows MA^', p even. 



4.3 Phase coexistence in the massive regime 

For the massive flow, a very different picture is expected at large distances. 
The ultraviolet perturbing operator is self-dual, and so in the massive direc- 
tion it should move the model onto a surface of first-order transitions. In 
this region the ordered and disordered phases coexist, and so the possible 
vacua should form a ^Ar-multiplet of ordered ground states, {a) = e~^, 
and a disordered ^Tv-singlet with (a) = 0. 

Finite-size scaling at a first-order transition is reasonably well-understood. 
In a cylindrical geometry, dominant configurations consist of a sequence 
of domain walls stretching across the 'spacelike' dimensions of the system. 
These lift the ground-state degeneracy, replacing it with an energy splitting 



which in d+1 dimensions is of order W^^"^^^ exp{—aR^) [45|. Here R'^ is 
the spatial volume and a the surface tension of the domain wall; in 1-|-1 
dimensions a is equal to the mass of the corresponding single kink. 

When there are m > 2 different infinite-volume ground states, the stan- 
dard discussion must be generalised a little. Consider first the situation 
where all kinks have the same mass, the vacua that they connect being en- 
coded in an incidence matrix lab- Then the usual instanton-gas argument 



(see for example ref. |Q) must be supplemented by a diagonalisation of I. 
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The prefactors of =Fl ~ minus the eigenvalues of o) ^ which distinguish 
the energies of the lowest-energy symmetric and antisymmetric states in the 
Ising model are replaced by — Aq < — Ai < . . . < — Am-i, the negatives of 
the m eigenvalues of /. (This can equivalently be understood via the effec- 
tive quantum-mechanical problem that remains once all transverse degrees 
of freedom have been integrated out.) In particular, the negative of the 
(unique) largest eigenvalue controls the leading behaviour of the ground- 
state energy [^]. Note however that all of the eigenvalues will be seen in 
the full spectrum, the others appearing in the asymptotics of those energy 
levels which become degenerate with the ground state in infinite volume. 
This explains the observation of one instance of this phenomenon made in 
ref. ||4^, and it offers the interesting prospect that it might be possible to 
'hear' the shape of the kink structure in the large- asymptotics of the de- 
generating energy levels. We checked that the excited-state TBA systems 



reported in refs. [48, ^ are consistent with this idea. 

The above assumed that all the kinks had the same mass. In general 
there will be n>l different kinds of kinks, with masses that can be ordered 
as M^(i)<M^(2)< • • • <M^(„) for some permutation vr of 1 . . . n. (In the W.n 
models, 7r(l) = n for A^<16, and 7r(l) = 1 thereafter.) For each mass Mi, an 
incidence matrix I^^ can be defined to encode the number of kinks of that 
mass which join each pair of vacua. If these matrices commute, then they 
can be simultaneously diagonalised to give sets of simultaneous eigenvalues 

,.(1) .(2) \(")^ > > Cx(l) \(") ^ 

[Aq ,Ao , • • • Aq ) > [A^ ,Ai ,...Ai ) ^ ■ ■ ■ ^ (A^n-i, Aj^_i, . . . A^^i) . 

(4.34) 

The ordering implied here is lexicographic, first by the eigenvalues Xp'^^^^ 
for the lowest kink mass, then by those for the second-lowest mass, and so 
on. The asymptotic behaviour of the k^^ energy level is then 

Ek{R) ~ -Al'^Ai(i?) - x'i'>A2{R) - ... - xPMR) , A: = 0, 1, . . . m-1 

(4.35) 

where Aj(i?), a function of order R exp(— Mji?), gives the leading energy 
splitting for a system with just two vacua and kink mass M,. This formula 
must be interpreted with some caution: if some of the heavier kinks have a 
mass more than twice that of the lightest, then their leading contributions, 
included above, will be less important than sub-leading contributions from 
the lightest kink, which have been omitted. Nevertheless, the functional 
forms of the Aj(i?) provide clear fingerprints for leading one-kink terms of 
each type, and match with the terms which emerge at the first step of an 
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iterative solution of the TBA equations. Thus while not always the leading 
asymptotic, ( [4.35 ) isolates exactly those terms that will be important below. 



The result seems to be consistent with a previously-known case. In 
ref. 1^^, excited-state TBA systems were proposed for the perturbation of 
the 2jv-symmetric conformal field theory by its first thermal operator e^^^. 
In contrast to the e^^^ perturbation that has been the principal concern of 
this paper, e*-^-* is anti-self-dual and moves the model off the self-dual hy- 
perplane. In the low-temperature direction, this leads to a theory of kink 
scattering. For current purposes it is convenient to group the kinks and an- 
tikinks together, even though their scattering is in fact diagonal, so that the 
model exhibits [N/2] different kink types, with masses Mj = M sin^nj /N), 
j = 1... [A^/2]. A number of different excited-state energies were found 
in ref. for those with conformal dimensions k{N—k)/2N{N+2) at 

short distances, we extracted the large R asymptotics in the way shortly 
to be described for the self-dual flows. The resulting eigenvalues were 
X^^^=2cos{27rjk/N), excepting the self-conjugate kink for even, 

for which A^^^^^= cos(7rA;). If it is assumed that these describe the excited 
states for untwisted boundary conditions in the low-temperature phase, 
then a compatible set of incidence matrices 1^^^ can be found on setting 
^ab ~ [-^N-ijab, the incidence matrix of the affine ajv-i Dynkin diagram, 
and then defining ij^ . . . ijj'^^'^^^ via an S'[/(2)-type fusion hierarchy: 

= jU-i) + /(i+i) , (4.36) 

with set equal to 21, twice the identity matrix. The factor 2 can be 
traced to the grouping of kink with antikink, which as before leads to an 
exceptional case for A^ even: the matrix /(^Z^) obtained from ( 4.36| ) must 



be halved. All of this is in agreement with the usual picture of the kink 
structure in the ordered phase, and in particular matches the bosonic parts 
of the supersymmetric solitons studied in [pO|]. 

There is a potential problem in more general cases: the assumption 
that the incidence matrices commute is rather strong, and one can envisage 
theories of kink scattering for which it fails. (For the example just described, 
it followed from the symmetry of the A^ degenerate vacua.) In the related 
context of integrable lattice models, the property is often built in from the 
start (see for example ref. [^]). However, from the S-matrix perspective 
it turns out to be a simple consequence of integrability, as shown by the 
following argument. 
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Consider the model on the full line, with spatial coordinate x, and restrict 
attention to the sector of those states which interpolate between the vacua a 
at X = — oo and c at x = +00. Irrespective of integrability, the conservation 
of topological charge implies that this sector is preserved under time evolu- 
tion. But if the model is integrable, then the sector can be further restricted 
to the (possibly empty) subspace of two-kink states with masses Mj 7^ Mj, 
and rapidities 9i 7^ 9j. If we take 9i > 6j, then for the in states, found as 
t — > — 00, the kink of mass Mj lies to the left of the kink of mass Mj. Hence 
Din, tti6 dimension of the space of in states, is equal to J2b^ab^bc- Now 
since Mj 7^ Mj and the model is integrable, reflection is ruled out and the 
space of out states is spanned by states in which the kink of mass Mj is to 
the right of the kink of mass Mj. Thus, Dout = J2b-^ab -^bc ■ since we 
can map from one space to the other and back by a unitary matrix S and its 
inverse, the two dimensions must be equal. That is. Din = Dout, which when 
taken for all pairs of vacua a, b is exactly the statement I^^^] = 0. Note 
that the argument holds even when some kinks have multiplicities higher 
than one, and also when there are some kinks of zero topological charge. 
Both situations will be encountered below. 

In this paper attention is being restricted to the ground-state energy. 
Even so, the formula (4.35) turns out to contain a large amount of informa- 
tion about the kink structure. To extract the asymptotics from the massive 
TBA equations, an iterative approach can be used [|l5|. To first approxima- 
tion the massive pseudoenergies behave as £f\0) ~ i^i^\9) = MiRcosh9, 
while the others take the constant values logXj^'^\ as given in appendix C. 
After one iteration, 

e[^\9) = M^Rcosh{9) + ^ N^f log [l + ^~^) , (4.37) 

where 

< = ^/^^<(^)d^ 

with ^'^^ a generic notation for the kernel linking e^f^ to in the TBA 



system. Hence 

E{R) = ^ / d^Mj cosh ^ log 

2'7r J— CO 



1-J2pd9 Mi cosh 9 e-A^'«-°^i^^ H + 



1/3 

(13)- " '■ 
j 



1 
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Identifying the factors Ai{R) = ^ d9 Mi cosh 9 e' ^^^^ , the eigen- 



values Aq ^ can be read off: 



n 

i,/3 



1/3 



(4.38) 



For the calculations it is actually easier to obtain the correction term in ( 4.37 ) 
by substitution into the relevant Y-system. The final results are remarkably 
simple: there are n different kink masses Mi, i = 1 . . .n, equal to the corre- 
sponding sine-Gordon masses ( |A.3| ), ( |A.4 ), and for the ground state energy 
we have 



A, 



A, 



(n) 



i+1 



1 



.n-1) 



(4.39) 



It remains to identify a set of commuting, non-negative integer val- 
ued, (A^-|-l) X (A^-|-l), symmetric and ^Ar-symmetric matrices I^^^ . . . I*-"'^ for 



which these form the largest (in the sense of equation ( 4.34 ) ) set of simulta- 
neous eigenvalues. In fact, the eigenvalues are individually maximal. To see 
this, first recombine the incidence matrices as Iab{R) = J2i ^i{^)Iab ■ ^ 
large enough, the eigenvector of Iab{R) with the largest eigenvalue is also the 
simultaneous eigenvector that we are looking for. Furthermore, for R finite 
Iab{R) is connected (note, the same cannot be assumed for the individual 
/^^■'). Hence, we can apply the Perron- Frobenius theorem to deduce that 
this eigenvector is ^^v-symmetric, with all components of like sign. Using 
the theorem in reverse, this eigenvector must be maximal for the individ- 
ual incidence matrices. Two of these matrices can be written down almost 
immediately. Up to a relabelling of the nodes 1 . . . corresponding to the 
ordered vacua, they are: 



I 


\ 




( 


1 \ 






An) _ 





1 


V ■ • 


2 J 




[ 1 1 . 


V 



(4.40) 



where ^at-i is the NxN incidence matrix of the affine a^^i Dynkin di- 
agram. So long as is not a perfect square, is fixed uniquely, while 
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for /(^^ we have to add a connectivity assumption when is not prime. In 
fact, any of the other possibihties for I^^^ would have the kinks of mass Mi 
joining non-adjacent ordered vacua. At least in a classical picture of the 
kinks, this can be ruled out since Mi is minimal among the remaining kink 
masses. There only remains the possibility for one of the blocks making up 
/^^) to be equal to zero, and imposing = rules this out as well. 

For the remaining kinks it is not possible to be so definite on the basis 
of the ground-state energy alone. However, one set of conjectures seems 
natural. Note first that the forms of I^^^ and /^"^ are consistent with the 
idea that the kinks of mass Mi are the (lowest-mass) bound states of a set 
of fundamental kinks with mass Mn=M. This parallels the related sine- 
Gordon model, where M„ is the mass of the soliton, and Mi the mass of 
its first bound state - the first breather. If the analogy is to extend, then 
the kinks of mass Mj in the model should be bound states of j of the 
kinks of mass Mi. The allowed topological charges should therefore be found 
among the paths with j steps on the graph of I^^^ . Adding this to the other 
requirements leads to the proposal 








. . . 


j+1 



(J 



l...n-ll 



(4.41) 



An alter- 



V 

Commutativity with /("^ was used to settle the value of l]P^i ^_,_i 
native characterisation is the statement that the breather-related incidence 
matrices form part of an 5'[/(2)-type fusion hierarchy: defining I^^^ = 1, the 
identity matrix, and I^^^ as in equation ( f4.40| ), the matrices I*-^-' . . . are 
determined by the relation ( [4.36 ), just as for the e^^^ perturbation, though 
with one more vacuum, and different initial conditions. Note that if I^^^ 
and commute and have the correct Perron-Frobenius eigenvalues, then 
the same automatically holds for /^^^ . . . given the fusion property. 

Some preliminary investigations of excited-state spectra for the ^n-related 
cases are also consistent with ([4.41| ). 

Figures § and ^ illustrate the general idea with the kink structures for 
A^=9 and A^=15. Tunneling to and from the disordered phase only occurs 
via the fundamental- kink instantons, while instantons associated with an 
increasingly intricate structure of higher kinks serve to connect the various 
ordered phases. Note also the presence of tadpoles in all but the fundamental 
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Figure 2: The matrices I^^^ and /^^^ for the massive theory 




Figure 3: The matrices I^^-* and I^^^ for the massive theory 



set of kinks. These have been seen before [52|, and their appearance in the 
finite-size effects can be understood either semi-classically, as corrections 
due to instantons of zero topological charge, or else as being due to virtual 
particle-like excitations ('breathers') above the various vacua. Either way, 
the predicted effect on the large- i? asymptotics is the same, and has indeed 



been observed in some other situations |47]. A more novel feature in this 
context is that, for A'^>16, some kinks appear in more than one incidence 
matrix from the set I^^^ . . . I^"^ . In other words, some pairs of vacua are 
joined not only by a simple kink of minimal mass, but also by excitations 
of this kink with higher masses - these can be thought of as breathers with 
non-zero topological charge. 

To see this aspect most clearly, the spectrum of asymptotic one-particle 
states can be reorganised according to their topological charges. Let a, a+i 
always label ordered vacua 1 . . . A^, with A^-|-l the disordered vacuum. With 
n as ever equal to the integer part of (A^— 1)/4, the prediction is that the 
massive flow should exhibit: 

• for i = 1,2,... n—1 : a kink joining each a to a+i with mass Mj , and 
excitations of this kink with masses Mj , j = i+2, 1+4, . . . < n—1 ; 

• kinks joining A^-|-l to each a with mass M„ ; 
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• excitations of each ordered vacuum with masses Mj, j = 2, 4, . . . < n— 1 ; 

• excitations of the disordered vacuum with masses Mi, M2, ■ ■ ■ M„_i . 

All appear with unit multiplicity at each mass except for those in the last 
line, which have multiplicities 2, 3, ... n. 

The discussion in this subsection has been based on a physical input to 
determine how the 2 tv symmetry acts on the vacua: we required one singlet 
and one A^-element multiplet. Whenever N factorises, say as N=PQ, then 



it is possible to use an orbifold construction [53| to produce alternative 
patterns of vacua, which might also be described by the self-dual TBA 
systems. The incidence matrices now split into blocks of sizes P and Q, 
rather than the previous N and 1. The basic matrices are: 





\ 




( 




[Ap-l\ab 





























it 






where ^p-i and ^q-i are the incidence matrices of the afhne apl^ and 
ciQ_i Dynkin diagrams, and lad is a PxQ matrix with all entries equal 
to 1. The remaining matrices, 

j(2) _ _ _ follow from 1^ 

via the fusion 

relation ( f4.36| ), with I^^^ = 1. Such systems of vacua and kinks could be 
relevant to first-order transition regions where 2p and 2q ordered phases 
coexist. 



5 Conclusions 

The checks performed in this paper have been rather exhaustive, and seem 
to us to establish beyond reasonable doubt that each pair of the self-dual 
2jv TBA systems listed in section 2 does indeed describe a pair of flows from 
the Zjv-symmetric conformal field theory into massless and massive phases, 
with the two flows being related by a change in the sign of a coupling con- 
stant. Given the earlier results of ref. |11], this effectively settles the question 



posed in the introduction about the location of the Fateev-Zamolodchikov 
Boltzman weights in the larger phase diagram. In addition, a number of 
predictions have been obtained which deserve further exploration. 
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One important outstanding question is to find a collection of S-matrices, 
both massive and massless, lying behind the TBA systems that have been 
proposed. In fact, a complete set of S-matrices for the massive flows was 
suggested some years ago by Fateev |17|, but we have not been able to 
reconcile these with the TBA results reported in section 4.3. In particular, 
Fateev predicted that for N odd there should be (N—5)/2 bound states 
of the fundamental kinks, all with different masses, while the asymptotics 
of the TBA systems proposed above predict only [(A^— 5)/4]. If the TBA 
systems are correct, it is very hard to see how about half of the kink masses 
could be lost from the instanton gas calculation. Given the stringent tests 
to which these systems have been subject in this paper, it is possible that 
the resolution will be found in some modification of Fateev's proposals. 

On a more mathematical note, the Y-systems and dilogarithm sum rules 
found here for odd appear to be new. Past experience (see for exam- 
ple ref. [54 1) would suggest the existence of new fermionic representations 
for the characters of the ^Tv-symmetric conformal field theories, and possi- 
bly also an alternative integrable deformation of the Fateev-Zamolodchikov 
Boltzman weights |8| which, contrary to that discussed by Jimbo et al |11], 
leaves the symmetry unbroken. 

Further study of the massless flows should be rewarding: as highlighted 
in section 4.2, the systems introduced in this paper add to an already-rich 
collection of flows arriving at c=l. The precise relationship between these 
models needs to be elucidated, and the possibility to extract at least some 
infrared information exactly should help in this task. Since the various 
infrared limits can all be regarded as irrelevant perturbations of a single free 
boson, they should also provide the simplest possible playground to explore 
more general issues in non-renormalisable field theory. 

Finally, there remain many unresolved features of the larger phase dia- 
grams of the 2 jv spin systems Q , beyond the particular self-dual directions 
out of the Fateev-Zamolodchikov points that have occupied our attention in 
this paper. We hope that these will provide further opportunities to apply 
the ideas of perturbed conformal field theory and the TBA to problems of 
independent statistical-mechanical interest. 
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P = r = ^r^- (A.2) 



A Sine-Gordon and kernel data 

The sine-Gordon model can be defined by the Lagrangian 

C = -{d^,ipf + ^cosP^ (A.l) 

and results in a family of scattering theories, conveniently parametrised by 
either p or h, where 

2 _ 

It will also be useful to define N = 2/i+4 = 327r//3^, and to set n equal to the 
integer part of (A'^— 1)/4. The spectrum contains a soliton-antisoliton dou- 
blet {s,s) and, for h > 2, n—1 soliton-antisoliton bound states, or breathers. 
If s and s have mass M, then the breather masses are 

Mk = 2M sin — = 2M sin , = 1, 2, . . . n-1 . (A.3) 

The soliton labels s and s will often be replaced by n and n+1, and accord- 
ingly 

Mn = Mn+l = M . (A.4) 

Scattering among solitons and antisolitons is generally non-diagonal; the 
amplitudes can be found in |55|. However, once a breather is involved the 
scattering becomes diagonal. By analogy with a notation used for the affine 
Toda theories, define the blocks 



sinh(f + ^) 



(^)W= ■ C/e ; {x} = (x-l)(x + l) . 

smhf'' i7rx\ 

Then 

Sjk = f{ {l}{h -1} (j, = 1 . . . n-1) , (A.5) 



li-fc|+i 

step 2 
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and 

h/2 + k-l 

Skn = Sk,n+i = i-l)'' n W (fc = l...n-l). (A.6) 

/i/2-fc + l 
step 2 

The couplings relevant for the 6„ and dn+i series of TBA systems cor- 
respond to h = 2n—l and h = 2n respectively - the dual Coxeter numbers 
of bn and dn+i- When h = 2n, the entire S-matrix is diagonal with the 
remaining S-matrix elements given by 



n even 



2n-3 


2n-l 

Sn,n+1 = Y\. ) 


step 4 


step 4 


2n-l 


271-3 


n w 




! = 1 
step 4 


;=3 

step 4 



In this case, the various prefactors of —1 (which anyway make no difference 
to the definitions below) can be omitted to leave the minimal version of the 
dn+i Toda S-matrix. 

The TBA systems in the text involve certain kernels. The first of these 
are related to the logarithmic derivatives of the diagonal sine-Gordon S- 
matrix elements: 

(pjk = logSjfc ; V'jfc = logTjfc . (A. 8) 

Here at least one of the indices j and k must have a 'breather' value, 
1,2, ...n—1. In the definition of ipjk, the function Tjk is obtained by re- 



placing each block {x} in (A.5-A.6) by (x). The kernels ■0jfc do not appear 
in the TBA equations for the sine-Gordon model itself, but rather they are 
part of the conjectures that are the main topic of this paper. Their forms 
are a natural generalisation of previously- known examples |]l^] . 

For the (i„+i-related models the definition can be extended immediately 
to cover the remaining cases when both j and k take the values n or n+1. 
Otherwise, the non-diagonal scattering of the solitons means that the asso- 
ciated kernels are more elaborate. Define an integer p hy N = An+p, and 
then set 

p cosh —6* 
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This function has the property that 



xp*f{0+'-^)+xp*f{o-'ii) = m. 

The sohton-sohton S-matrix consists of a 'scalar' piece multiplying a non- 
diagonal matrix. Associated with the first factor is a pseudoenergy SniO) , 
with interaction kernels in the TBA given by 

4>nn(.0) = Xp*4>n,n-l{0) , i^nn{0) = Xp*i^n,n~l{0) ■ (A.IO) 

(For N<7, (j)n,n-i and V'n.n-i are not defined and we set (t>nn='4'nn=^-) In 
addition, depending on the continued- fraction expansion of 2//i, the non- 



diagonal part induces a number of fictitious 'magnonic' particles |56], to- 



gether with their attendant kernels. We will only need those which arise 



when N is an integer; these can be obtained from the Y-systems of ref. |24], 
and are 



cosh h9 



0) 



cosh 2h9 



cosh 



8/1 cosh 
^'^^^ - 3(4cosh2f^-3) 

8/icoshf0 , 



As above, h is equal to N/2 — 2 . 



B Y-systems 

This section records the self-dual 2jv Y-systems. Their derivation from the 
TBA equations is rather involved, but follows the same lines as the simpler 



cases explained in ref. [21|. In particular, the energy terms are eliminated 
using the following variant of the Perron-Probenius mass property of the 
ade-related theories: 

2cos(f)Mi = (i = l,...n-2); 

i 

2cos(f)M„_i = M„_2 + 2cos(i^)M„; 
2cos(|f)M„ = M„_i, 
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where = An+p, as in the last section, and is the incidence matrix 

of the a„-i Dynkin diagram. 

The results are 22-symmetrised versions of the sine-Gordon Y-systems 
constructed in ref. [^4|. Figure ^ shows a diagrammatic representation, 
following refs. |18, 24|. In order to write them down in a managable form, 
define 

Y\"\e,r] = y/°)(0-lH:)yW(0 + Y) 
Yf'\9,r} = (i + y/")(^-if:))(i + y/")(^ + Mi; 

With this notation in place, the Y-systems for the four cases can be written 
as follows. 

1) iV = 4n+l 

Nodes {i<n—l,a = 1,6), with a = 7— a : 



Nodes (n— l,a=l,6) : 
Y^^\[9,l] = (l + Y^ 



i%{er'y\i + Y^i'Me)) (i + Y('\9)) (i + yi^)(^)) 



xyi^){e,3/4}yf {^,2/4}yf {^,1/4} 



Y^:UoA] = {i + Y^'Me)-Y{' + ^n%io)){i + Y(^\9)){i + Y('He) 

Nodes (n, a=l ... 6) : 



X yf {^, 3/4}yf {^, 2/4}yf {^, i/4} 



eel 



with yi"^(6l) = when a / 1, 6. 



2) N = 4n+2 (bn) 

Nodes (i<n— 1, a=l, 3) : equation (|B.1|) holds with q = 4— q . 
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1 2 n-1 n,l 




Figure 4: Diagrammatic representations of the Y-systems. The 

left/right columns represent the massive/massless theories respectively. 
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Nodes (n— l,a=l,3): 

Nodes (n, a=0 ... 4) (with ¥^"^{{9) = when a / 1, 3) : 

Y^:\eM2] = {i + Yt\{e))\{{i + YP{e)-Y"' 

3) iV = 4n+3 

Nodes (i<n— 1, a=l, 3) : equation ( |B.1| ) holds with a = 7— a . 
Nodes (n— 1, a=l, 6) : 

y^nUoA] = {i + Y^'Me)){i + Y('\6)){i + Y^%{e)-'y'Y^^\e,i/A} 
y^n\[o,i] = {i + Y^%ie)){i + Y^'\e)){i + Y^'Mo)-'y'Yi':\e,i/A} 

Nodes (n, a = 1 . . . 6) : 

Y^^\e,3/A] = {^i + Y}^\{e)){i + Y('He)-'y\i + Y('\e)-'y' 

xV('He,2/A}vP{9,l/4} 

Ff[e,i/4] = {i+Y:i'\0)){i+Y('\e)-'y' 

Ff [0,1/4] = (i + yj^) {6)){i + yi^) (0)) (i + yjs) (e)) 
[0,1/4] = (i + yp)(0)) 

F(f)[0, 1/4] = {^i + YP{e)){l + Y(^\9)-'y' 

Ff[0,3/4] = (i + ySw)(i + yi^)(0)-^)"'(i + yi^)W-^)"' 

xVP{9,2/A}vP{e,l/4} 

4) N = 4n+4 

In this case the scope of equation ( [B.lD , with 5 = 3— a, can be extended to 
cover all of the nodes {i=l . . . n+1, a=l, 2) : 

fS")[0,i] = (i + y/")(0ri)" n(i + ^i"^w)" 
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In the above formulae, ?|^" ^\ /|^|^ and are the incidence matrices 

of the corresponding Dynkin diagrams. 

We checked numerically the periodicity implied by these Y-systems up to 
A^=30, and found 

which in turn implies the result (3.7) quoted in the main text. 



C Dilogarithm sum rules 

The sum rules are conveniently given in terms of certain limits of the func- 
tions Y^"\6) = e^i With the £f^\9) the solution of a massless system, 
define 

T^"^ = lim K^^Ve) ie finite) , 

* Affile ^ ^ ' 

x!>"^ = lim y/"^(6') (ikffi finite), 

= hm Yt\e) (0 finite). (C.l) 

MR^ 00 

(For the pseudoenergies which solve the massive systems, the T^"^ and X^'^'^ 



are unchanged while 2^°^^ = X^"^^ .) The numbers defined by (C.l) furnish 



stationary solutions to the Y-systems of appendix B, subject to the following 
constraints: all T^"^ are finite; X^^^ = 00 Vi with the rest finite; and z\^'^ = 

= 00 yi with the rest finite. The infinite quantities do not contribute 
to the sum rules. To find the values of the others, start with the following 
ansatz for the On-i ^ cl2 t^-il^ 

r^{a) _ sin((i + 3)7]) sm{irj) 
* sin((2?7)) sin(ry) ' 

Xp = {i + 2)i, 

with i = 1 . . . n—1. The value of rj can be found by examining the final nodes 
(n, a); a positive solution results on setting 77 = jy+l ■ This also gives values 
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in good agreement with the direct numerical solution of the TBA equations. 
For the remaining nodes the expressions are more complicated, and will be 
given case-by-case. 

1) = An+1 



sin((2n + 2)?7)sin(2??) ' " " 1 + T^^l^ - {r'h^ 

^(4) ^ (Ti^))3/2(i^^(l))l/2 

^(3) ^ (t(^)3(1 + t(^)) 

1 + 2Ti'2i + (Ti'2i)2 _ 2(t(,'))2 - 2Ti'2i(T(,'))2 - (t(/))3 



2n + 1 ' " (3n + l)(n-M) 
= A'?) = , ;f(3) 



2) iV = 4n+2 



3n + 1 ' " 8n2 + 1 + 6n 

^(4) _ ^(2) _ 7(5) - i ^(3) _ i 



sin2(n?7) sin(n?7) 



sin((2n-F2)?7)sin(2ry) ' " " sin((n 2)r/) 
T(2) = 



sin^(n77) 



sin^((n + 2)??) 



3) iV = 4n+3 



2n-M ' " " n + 1 ' " 



sin^ (nr/) 



^(1) = T^^) = 

" " sin((2n + 2)r/)sin(27?) 

sin((n + 2)77) sin((n + 1)77) + sin2((2n + 2)r/) 
sin2((n + 2)r/) 

-^(3) _ ex (^))2 _ n ^(2) _ -f-(5) _ rt^)^ - 1 _ , 
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v(6) ^ ;p{5) ^ n{3n + 2) 

2n + 1 ' " (n + 1)2 
(3) ^ (4n + 3)(2n + l) .4) ^ (2) ^ 3n + 2 
(n + l)2 ' " n + 1 

4) N = 4n+4 

T(") = t(") = , , , {a = 1,2). 

2sin(r?)cos((n+l)??) ^ ' ^ ' 



n 



These numbers enter into the central charge calculation via the arguments 
of Rogers' dilogarithm functions 



1 



Iny _^ In(l-y) 



.1-2/ y 

The following sum rules were verified numerically, for all cases up to A^=30 : 

4n+l 



Cx 




2(jV-l 

N + 2 ' ] 2 for iV = 4n+3 

4n+4 



cx 









for N = 


-\ 


5/2 


for iV = 






for N = 




[? 


for N = 


4 


for 


= 4n+l 


5/2 


for 


= 4n+2 


2 


for N 


= 4n+3 


1 


for N 


= 4n+4 


3 


for iV 


= 4n+l 


3/2 


for N 


= 4n+2 


1 


for iV 


= 4n+3 





for iV 


= 4n+4 



In fact, the identities found were more general than this: in the spirit of 
ref. [^], we checked that they also hold if the constants X, T and Z are 
replaced by any non-stationary solutions to the Y-systems of the last section, 
so long as an additional average over a full period is performed for each 
result. These periods are 2(A^+2) for the T's, for the A:"s, and 8, 4 and 
8 (when p = 1,2 and 3) for the Z's. 
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The sum rules are not obviously relevant when N=An+3, since for these 
TBA systems the kernels given in the text are not all symmetric. However, 
at the expense of some additional complexity, the systems can be rewritten 
in a symmetrical way, after which the central charge calculation goes through 
as usual. The forms given earlier were selected as the simplest versions which 
are also suitable for direct numerical solution. 

With this detail out of the way, the ultraviolet and infrared central 
charges are 

CUV = CT-cx= (^pj^2) ' ^ 

CIR = cx - cz = 1 . (C.3) 



D Regular UV expansion coefficients 

This appendix tabulates fits of the numerical solutions to the TBA equations 
at small R to the expansions 

oo 

RE{R) + (bulk piece) = 27r V {MRf^"" . 

6 n 

The scale M was set equal to 1, and c to 2(A^— l)/(A^+2). The bulk pieces, 
and some details of the numerical procedures used, can be found in section 3 
of the main text. For completeness, coefficients for N=8, obtained simply by 
doubling those reported by Zamolodchikov for the tricritical Ising to Ising 
flow in ref. have also been included. 
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N=5: 



m 


Fm (massive) 


Fjn (massless) 





1.98 X lO"^'^ 


-3.86 X 10"-'-' 


1 


-7.78 X 10-^2 


3.23 X 10-^2 


2 


0.0112031704 


0.0112031702 


3 


-3.37 X 10-°9 


1.35 X 10-°9 


4 


0.00038484 


0.00038477 


5 


0.019215 


-0.019217 


6 


0.02056 


0.02055 


7 


0.0001 


-0.0006 



A^=6: 



m 


Fm (massive) 


Fm (massless) 





2.15 X 10-18 


4.10 X 10-1^ 


1 


1.01 X 10-13 


1.22 X 10-13 


2 


0.026663544239 


0.026663544230 


3 


0.0202428666 


-0.0202428658 


4 


0.00088616 


0.00088614 


5 


-0.0061087 


0.0061090 


6 


-0.002800 


-0.002803 


7 


0.00033 


-0.00031 
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N=7: 



m 


Fm (massive) 


Fjn (massless) 





3.98 X 10-^5 


8.02 X 10-1^ 


1 


2.42 X 10-1'' 


-5.84 X 10-15 


2 


0.04302032556 


0.04302032557 


3 


0.0128965993 


-0.0128965995 


4 


0.000972018 


0.000972029 


5 


0.00120282 


-0.00120286 


6 


0.0000682 


0.0000690 


7 


-0.0000085 


0.0000077 



A^=8: 



m 


Fm (massive) 


Fjn (massless) 


2 


0.06589308648 


0.06589308648 


3 


-0.02744520208 


-0.02744520208 


4 


0.00073623262 


0.00073623264 


5 


0.00102222 


-0.00102216 


6 


-0.00016842 


-0.00016851 


7 


-0.00004770 


0.00004710 


8 


-0.0000268 


-0.0000226 
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A^=9: 



m 


Fm (massive) 


Fm (massless) 





7.27 X 10-^-^ 


2.89 X 10-^3 


1 


-4.21 X 10-^2 


-1.97 X 10-^1 


2 


0.124757776 


0.124757776 


3 


-0.01159266 


0.01159265 


4 


0.0004886 


0.0004887 


5 


0.000574 


-0.000575 


6 


0.000088 


0.000089 


7 


-0.000028 


0.000025 



N= 10: 



m 


Fm (massive) 


Fm (massless) 





1.17 X 10"^^ 


-2.94 X 10-1^ 


1 


-8.58 X 10-^^ 


2.40 X 10-^2 


2 


0.0247185731 


0.0247185730 


3 


-0.006211643 


0.006211644 


4 


0.00031116 


0.00031115 


5 


0.00030406 


-0.00030402 


6 


0.0000402 


0.0000401 


7 


-0.0000137 


0.0000138 
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